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ABSTRACT 

This paper addresses the challenge of understanding the typical star formation histories of red 
sequence galaxies, using linestrength indices and mass-to-light ratios as complementary constraints 
on their stellar age distribution. We first construct simple parametric models of the star formation 
history that bracket a range of scenarios, and fit these models to the linestrength indices of low- 
redshift cluster red-sequence galaxies. For giant galaxies, we confirm the downsizing trend, i.e. the 
stellar populations are younger, on average, for lower a galaxies. We find, however, that this trend 
flattens or reverses at a < 70km s~^. We then compare predicted stellar mass-to-light ratios with 
dynamical mass-to-light ratios derived from the Fundamental Plane, or by the SAURON group. For 
galaxies with ct ~ 70 kms"'^, models with a late "frosting" of young stars and models with exponential 
star formation histories have stellar mass-to-light ratios that are larger than observed dynamical mass- 
to-light ratios by factors of 1.7 and 1.4, respectively, and so are rejected. The single stellar population 
(SSP) model is consistent with the Fundamental Plane, and requires a modest amount of dark matter 
(between 20% to 30%) to account for the difference between stellar and dynamical mass-to-light 
ratios. A model in which star formation was "quenched" at intermediate ages is also consistent with 
the observations, although in this case less dark matter is required for low mass galaxies. We also find 
that the contribution of stellar populations to the "tilt" of the Fundamental Plane is highly dependent 
on the assumed star-formation history: for the SSP model, the tilt of the FP is driven primarily by 
stellar-population effects. For a quenched model, two-thirds of the tilt is due to stellar populations 
and only one third is due to dark matter or non-homology. 

Subject headings: galaxies: clusters: galaxies: elliptical and lenticular, cD: galaxies: stellar content: 
galaxies: fundamental parameters: dark matter: surveys 



1. INTRODUCTION 

Red-sequence galaxies (RSGs) dom inate the stellar 
mass in clusters and comprise ~ 70% (iBell et all 120031) 
of the stellar mass in the Universe. However, other 
than the fact that the bulk of their stellar population 
is older than ~ 1 Gyr, a detailed understanding of 
their star formation histories remains elusive. RSGs fol- 
low a tight color-magnitude r elation (hereafter C MR, 
ISandage fc Visvanathaiil 119781 : iBower et al] Il992f ) and 
scaling r elations such as the Fundamental Plane (here- 
after FP. iDressler et al.lll987l : iDjorgovski fc Davislll987f l 
with little scatter. At the same time, the dependence of 
the colors on magnitude suggests that there are indeed 
variations in their stellar properties as a function of lumi- 
nosity or a mass-related parameter such as stellar mass or 
velocity dispersion. Variations in the stellar populations 
would also affect the "tilt" of the Fundamental Plane 
with respect to the virial scaling, although the degree 
to which the tilt is due to stellar populations compared 
with dark matter (DM) and non-homology remain s hotly 
debated (t Prugniel fc Simien' '1996'; 'Pahre et al.' [199^ 
Gerhard et al. 2001; Truiillo ct al. 2004; CappcUari et all 



be b r oken by using stell ar absorption lines (jWorthevI 
Il994t iThomas et aLll2003f ). Recent results suggest that 
metallicity, a- enhancement and age vary alon g the 
mass or velocity dispersion sequence fC aid we ll et al.l 
[200l iNelan et al.l [20051 : iThomas et a l. 2005), an d also 
vary as a function of env ironment ( Thomas et al.l [20051 : 
iSmith et al.[[200l [2009bl) . 

A limitation of the above analyses is that, in gen- 
eral, the absorption line indices are sensitive only to 
weighted mean stellar age, with a weighti ng that strongly 
favor s the most recent star formation ^Serra fc Trageil 
|2007[) . The indices provide only weak information about 
the range of stellar ages within a galaxy. For conve- 
nience, then, stellar population ages are often quoted as 
the equivalent age of a Single Stellar Population (here- 
after SSP). In pr actice, other scenari os have been stud- 
ied. For example, iTrager et al.l (|2000D considered "frost- 
ing" models in which a small (few percent) "icing" of 
young stars on top of a "cake" of older stars yields 
a you ng luminosity-weig hted age; IThomas et al] (20021 
|2005[) modeled star f ormation histories as a Gaussian in 
age: iBell et al.l (I2003D assume an exp onential star forma- 
tion htstoryTKau^man^^Ein] (|2003l ) model SDSS galax- 
ies as a combination of an exponential with small bursts; 
and lHarker et al.l (|2006) considered models in which star 
formation starts early and is later quenched. In terms 
of predicting the line indices, these scenarios are largely. 



20061 iLa Barbera et al.ll2008l) . 



Broadband optical colors are not good discriminants 
of stellar p opulations beca use of the age-metallicity 
degeneracy ([WorthevI [l99l . These de generacies can 
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but not completely, degenerate, as we will show in Sec- 
tion [Hill 

However, it may be possible to break these degenera- 
cies using other, complementary, observations. One ap- 
proach is via chemical abundances: because different el- 
ements are produced by different progenitors at various 
times, it is poss ible (in principle) to use the abundances 
as a clock (e.g. iThomas et aLll2005[ ). In practice, how- 
ever, there are a large number of unknowns (masses of 
progenitors, yields, gas inflows and outflows) which ren- 
der the quantitative conversion of element abundances 
into time-scales problematic. An alternative approach is 
to compare the predictions of various star-formation his- 
tory scenarios against other observational relations be- 
sides line indices. At low redshifts, the observables in- 
clude the slope of the CMR, the tilt of the FP and the 
scaling of mass-to-light ratios with mass or velocity dis- 
persion. At higher redshift, one can compare the predic- 
tions for the evolution of the zero-point, slope and scatter 
of the CMR and FP, and of the RSG dwarf-to-giant ratio. 

In this paper we fit the observed linestrength in- 
dices adopting six generic model s of the star for- 
mation history Pr evious studies (jTrager et al.l l2000l : 
iBernardi et al.1 [20051) find velocity dispersi on to be the 
"drivin g parameter" of stellar populations. ISmith et al.l 
(I2009ar) . using the same spectroscopic data as in this pa- 
per, not only confirm velocity dispersion to be the pri- 
mary parameter of stellar populations, but find no ad- 
ditional dependence of the stellar populations on stellar 
mass. Therefore, for each model we assume a mean scal- 
ing and intrinsic scatter in metallicity, a-enhancement 
and "age", all assumed to be functions of velocity dis- 
persion. Here, "age" is some parameter related to the 
timescale of star formation, the definition differing from 
model to model. Having constrained the scalings via the 
line strength data, we generate synthetic clusters using 
the stellar population parameters to determine mass-to- 
light ratios, colors, magnitudes etc. For these synthetic 
clusters, we construct the CMR and FP relations and 
compare these to observational data from low redshift 
rich clusters, under the assump tion that these relations 
are universal for rich clusters ijLopez-Cruz et al] 120041 : 
iMcIntosh et al.ll2005f ). In a future paper, we will extend 
this approach to comparisons with high redshift clusters. 

An outline of this paper is as follows: Section [5] de- 
scribes the spectroscopic and photometric data sets; Sec- 
tion [3] describes the models of star formation history. 
The results of the fits to these models are presented in 
Section [4] Having fixed the parameters of the models 
with the spectroscopic data, we create synthetic clusters 
based on these models in Section [5] and compare these 
to dynamical mass-to-light ratios, colors, and the Faber- 
Jackson and Fundamental Plane relations in Section [6l 
We discuss the effect of systematics, the IMF and of dark 
matter on our results in Section [7| and discuss the impact 
of our results on derived stellar mass densities and the 
tilt of the FP in Section [H We summarize the results in 
Section [91 

Throughout this paper, we assume the following 
cosmological parameters: (f^M, f^A, h) = (0.3, 0.7, 0.7). 



2. DATA 



In this paper, we will make use of two datasets for 
the RSG population in rich clusters. We assume that 
RSG populations in rich clusters are universal, in the 
sense of having the same distribution of star-formation 
histories at a given RSG mass. Absorption linestrength 
data from deep observations of the RSGs in three rich 
clusters in the Shapley supercluster are used to derive 
the ages and metallicities of the RSG population. The 
resulting predicted colors, magnitudes and Fundamental 
Plane parameters for the RSG population are compared 
to data in the Coma cluster. 

2.1. Linestrength Data from Cluster Galaxies in 
Shapley Concentration 

The Shapley cluster sample con sists of ^ 180 R 
< 18 emission-free cluster galaxies ([Smith et al.l |2007|) 
that cover a wide range in central velocity dispersion 
(38kms~^ < a < 313kms~^). There is no morpho- 
logical selection in this dataset: the only selection is on 
Ha emission, and this is strongly correlated with color: 
~ 99% of the emission-free galaxies are red. Further- 
more, red-selection and morphology are tightly corre- 
lated: for example, only ^ 20% of NFPS galaxies with 
fj > 70 kms~^ are types Sa or later (Hudson et al. 09, 
in preparation). Finally, we note that nearly all of the 
galaxies in the Shapley cluster sample lie within half the 
virial radius (r2oo /2) of the center of one of the three rich 
clusters (A3556, A3558, A3562) in the supercluster. 

The Lick linestrength indices were measured from 8- 
hour, high signal-to- noise ratio (~ 60 A~^) spectra, 
obtained with the AAOmega instrument at the 3.9m 
Anglo-Australian Telescope. The spectra were obtained 
through a 1" radius fiber, corresponding to a phys- 
ical radius of 0.95 kpc. In this paper, we use six 
Lick linestrength indices (three Balmer indices and three 
metallic indices) to measure age-related quantities and 
metallicities, as discussed in Section [Jj below. 

2.2. Coma Cluster Calaxy Data 

We will compare predictions from our models to pho- 
tometric data fo r 218 galaxies in the C oma Cluster 
from SPSS PR 6 (IAdelman-McCarthvl l200l and 2MASS 
([Skrutskie et al.l l2006f ) with measured velocity disper- 
sions from the literature. We compare with Coma rather 
than Shapley itself because no published photometry is 
available for Shapley, whereas a large range of wavebands 
is available in Coma {u through K). We note that the 
Coma clu ster is of similar r i chnes s to the Shapley clus- 
ters, and iLopez-Cruz et all (|2004f ) have shown that for 
Abell clusters, the cluster-to-cluster dispersion in CMR 
colour (at the characteristic luminosity L*) is only 0.05 
\ii B — R. Moreover, the CMR of Coma is typical, deviat- 
ing by no more than 0.05 mags from the average relation. 
Figure [1] compares the line strength indices fr om the ho- 
mogeneous NFPS survey ([Nelan et al.ll2005f l. There is 
no significant difference between the index-cr relations of 
the two clusters. 

Derived photometric parameters in Coma include half- 
light radii, surface brightnesses within the half-light ra- 
dius, and total magnitudes measured in the SDSS r- 
band. These parameters have been corrected for the 
effects of seeing by fitting a model using GALFIT 
([Peng et al. 2002) and deriving corrections to observed 
values from the difference between the seeing-convolved 



Star Formation Histories of Red Galaxies 



3 



CO. 

X 




o 




-> ° 














,0 




















1.4 1.6 1.8 2.0 2.2 2.4 2.6 
log(o) 



2.8 




1.4 1.6 1.8 2.0 2.2 2.4 2.6 
log(o) 



2.8 




1.4 



1.6 



2.0 2.2 
log(o) 



2.4 



2.6 



2.8 



1.4 



1.6 



2.0 2.2 
log(o) 



2.4 



2.6 



2.8 



Fig. 1. — Linestrcngth-cr relations from the NFPS sample for Coma (blue) and the Shapley clusters (red). NFPS data for other clusters 
are shown in grey. In all cases, galaxies within 0.5 Mpc are used. There is no significant difference between Coma and Shapley clusters. 



model and the unconvolved model. In order to com- 
pare colors with ages and metallicities predicted from the 
Shapley linestrengths, we measure the colors within a 2" 
radius aperture corresponding to a physical aperture of 
0.97 kpc, which matches the 1" aperture of the AAOmega 
data in Shapley. (Shapley is approximately twice as dis- 
tant as Coma.) These aperture magnitudes are corrected 
for seeing in the same way as de scribed above . Dat a 
are k-corrected using formulae from lFrei fc Gunnl (I1994D, 
and co rrected for Galactic extinction using lSchlegel et al.1 
()1998l ). Finally, because there is evidence for gradients 
in the s tellar populations as a fu nction of cluster-centric 
radius (iSmith et al.ll2006ll2008bD . we limit fits to Coma 
galaxies within r2oo/2 (i.e. within 1.24 Mpc, 0.71°) of 
the cluster center. 

In order to model stellar masses, as well as to calcu- 
late aperture corrections, we will require the relationship 
between effective radius and velocity dispersion. This re- 
lationship in Coma is shown in Figure [51 We find that a 
broken power-law is required to fit the data: 
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logi?e = 2.715 logo- -5.81; 



log(cr) < 2.23 (1) 
log(CT) > 2.23 
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where i?e is in units of kpc and a is in units of km/s. 
The break in Re is closely related to the well-known 
break in the surface brightness behavior (jKormendvl 
1 19851), and has been seen in t he Faber- Jackson rela- 
tion (jMatkovic fc Guzmanll2065l ) and the Re — L relation 



Fig. 2. — Coma log(i?e) — log(a') relation. A broken power law 
(Equation |2]l was fit to Coma galaxies with cluster-centric radius 
less than r'20o/2 (filled circles) to match the Shapley sample. 
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iHvde fc Bernardil (|2008af) . The dispersion in log(i?e), is 
0.213 for log(cr) < 2.23, and is 0.251 for larger values. 
The scatter at the low-cr/faint end is comparable to that 
found at fixed luminosity for t he SPSS for fain t early- 
type galaxies; for comparison (jShen et al.ll2003l ) find a 
dispersion > 0.22 in log(i?e) for Mr > -19. 

3. SYNTHETIC LINESTRENGTHS FOR GALAXIES WITH 
COMPLEX STAR FORMATION HISTORIES 

Our goal is to extend the SSP analyses of previ- 
ous work to more complex star formation histories 
(CSFHs). A CSFH model can be generated by con- 
volving the SSP response with the star form ation rate. 
Here we co mbine the SSP SED's of Maraston ( Marastonl 
[1998', 2005) with t he a— enhanced Lick-in dex absorption 
linestrengths ( Thomas et al.l |2003| . |2004 together here- 
after the TMBK models) to construct linestrength in- 
dices for CSFHs. Specifically, the TMBK line indices are 
based on SSP models. To obtain predictions for complex 
star formation scenarios, we reconstruct the fluxes in the 
sidebands and in the central (absorption) band for each 
step of star formation. By integrating these fluxes over 
all time steps, appropriately weighted, we can produce a 
set of indices for any formation scenario. A very similar 
method applies to mass-to-light ratios. 

The range of possible star formation histories is clearly 
large and highly uncertain. Here we consider six sim- 
ple star formation scenarios that are intended to bracket 
more complicated, realistic scenarios. All but one of the 
models we consider have a single age-related quantity, 
whose definition depends on the specific model detailed 
below. The exception is the Old SSP model, which has 
a single fixed age. Note that "ages" from Lick indices 
have poor age sensitivity at old ages, so it is impossible 
to discriminate, for example, an old instantaneous burst 
from an old burst of the same age, but with a short du- 
ration of 1 Gyr. The models and their abbreviations are 
as follows: 

• Old SSP (hereafter OSP): a single burst of star 
formation at a fixed age of 13 Gyr. 

• Single Stellar Population (SSP): a single burst of 
star formation, with fitted age ^ssp- 

• Exponential star formation rate (EXP): star for- 
mation begins 13 Gyr ago, ends 0.1 Gyr ago, with 
a fitted exponential decay parameter r. 

• "Abruptly" Quenched (AQ) star formation: a 
model with constant star formation beginning 13 
Gyr ago, and an abrupt truncation of star forma- 
tion at a fitted lookback time, ^aq- 

• A 2-component "frosting" model (FR) consisting 
of a dominant (98% by mass) 13 Gyr SSP and a 
secondary (2% by mass) burst with fitted age ^fr- 

• A "Strangulation" model (STR): a quenched (con- 
stant SFR) model beginning 13 Gyr ago up to fitted 
tgxR, followed by an exponential SFR from <str to 
0.1 Gyr ago with fixed r = 1 Gyr. The mass frac- 
tions in each component are a function of tsTR, 
such that the SFR at isTR is continuous. 



In each model, we assume a single metallicity and a 
single a-enhancement for all generations of stars, and fit 
these parameters to the data. Our models therefore do 
not attempt to include chemical enrichment in any phys- 
ically motivated way and so the abundance parameters 
are intended to represent luminosity-weighted mean val- 
ues. 

Note that in the OSP, AQ, EXP, FR and STR models, 
we have chosen to fix the age of the old population or the 
"starting time" for star formation at 13 Gyr. While, in 
principle, one would like to leave the starting time as a 
free parameter, in practice, the multiple age parameters 
become highly degenerate. We will discuss the implica- 
tions of varying the age of the oldest stars in Sections 
[OJandO 

These models can be subdivided into two classes. In 
one class of models (OSP, EXP, FR), the stars are pre- 
dominantly old. In the case of the EXP and FR mod- 
els, intermediate ('^ 5 Gyr) luminosity-weighted ages can 
arise, but these are due to a small fraction (by mass) of 
very young (< 1 Gyr) stars in addition to the dominant 
(by mass) old population. In the other class of models 
(SSP, AQ, STR), there is little star formation at late 
times and the luminosity-weighted intermediate ages are 
due to stars which indeed formed at intermediate times. 
We will show that these two classes have different prop- 
erties, particularly in terms of their stellar mass-to-light 
ratios, even when constrained to match the same line 
index data. 

4. FITTING STAR FORMATION HISTORY MODELS TO 
LINESTRENGTH INDICES 

The goal of this paper is to create synthetic clusters of 
simulated galaxies with distributions of age and metal- 
licity that are consistent with the observed line indices. 
While the observational data are of high S/N, the fitted 
ages and metallicities of individual galaxies have sub- 
stantial correlated errors, particularly for the faintest, 
low velocity dispersion galaxies. Thus, rather than in- 
vert the data, we instead model the mean and scatter 
in the stellar population parameters, P ("age", metal- 
licity and a-enhancement) that are required to match 
the median linestrength and the scatter in the observed 
linestrengths. 

By comp aring the scatters in either linestrengths 
(|Trager et a l. 2000), colors (Bc rnardi et al . 2005) or to- 
tal M/L (|Cappellari et al.1l2006h as a function of veloc- 
ity dispersion, with the same scatters as a function of 
luminosity, one finds that the correlations with velocity 
dispersion are always tighter than those with luminosity. 
Indeed, we have confirmed this using our Coma cluster 
data described above. These results strongly suggest that 
the fundamental or "driving" pa rameter of stellar p opu- 
lations is the velocity dispersion ([Smith et al.ll2009al ). an 
assumption we will make in this paper. Thus we divide 
our sample into five bins by log(cr) and measure median 
Lick indices and scatters for each bin, both of which are 
then fit to each star formation history model. 

We will use the six Lick indices given in Table [1] to 
break the degeneracies between "age" (traced primar- 
ily by the Balmer lines: H(5F, H7F, H/3), Fe (Fe4383, 
Fe5015) and a-element enhancement (Mgb5177). We 
assign the galaxies to five bins in velocity dispersion, 
and for each bin we calculate the observed median 
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TABLE 1 
Best-Fit Line Index 
Corrections 



Line Index 


Correction" (A) 


H<5F 


-0.063 


H7F 


+0.026 


Fe4383 


+0.257 


H/3 


+0.061 


Fe5015 


+0.667 


Mgb5177 


+0.274 



" Applied to all observed bins 



linestrength index. Our stellar population parameters 
include one age-related parameter (tssp,tAQ,tsTR,T,tF-R), 
[Z/H] and a/Fe. 

In general, however, we expect that the linestrength 
data may have small zero-point offsets. To c o mpare with 
previous results, and following 'Smit h et al.l (|2009af l. we 
force the median line indices in the highest-velocity dis- 
persion bin to equ al the SSP-parameters recovered by 
iNelan et al.l (|2005f ) in their highest velocity dispersion 
bin: (age, [Z/H], [a/Fe]) = (10.8 Gyr, 0.24, 0.28). This 
corresponding offset given in Table [1] in each line is then 
applied to all bins. 

Our goal is not only to fit the median line indices at 
each velocity dispersion, but also fit the spread in ages, 
metallicities and a- enhancement to the spread in each 
line index of each bin. We first calculate the total ob- 
served scatter Stot in a given index / at a fixed cr by 
measuring the semi-interquartile-range (SIQR) and con- 
vert this to its Gaussian equivalent: S'tot = SIQR/0.67. 
This is more robust to outliers than the usual root-mean- 
square. We then estimate the intrinsic dispersion in line 
index properties as follows: 

where S'mcas is taken to be the median measurement error 
in the bin for each specific absorption line. 

4.1. Results 
4.1.1. Fits to Median Linestrengths 

We employ a simple fitting scheme to the observed 
median linestrength indices. The results are shown in 
Figure [3] and are tabulated in Table [21 along with the 
Shapley binned values to which the models were fit. Cor- 
responding model parameters appear in Table [S] Note 
that, for the Balmer lines H7F and 11/3, the scaling 
is roughly linear at high velocity dispersions, but the 
lowest-(T bin appears to deviate from this linear relation- 
ship. This will yield older ages for this bin. We will 
discuss this further in Section [5T4l 

From Fig. [H we see that all models, except the Old 
SSP Model which is clearly a poor fit, yield similar pre- 
dicted linestrengths. However, closer inspection of H/3 
shows the SSP, AQ and STR models are a better fit to 
this index. To quantify this, we calculate a for each 
model, for each line, as well as an aggregate over all 
six lines for the four lowest mass bins (since we have fit 
a free offset to match the high velocity dispersion bin), 
and these values are tabulated in Table [21 Formally, the 
EXP, FR and OSP are poor fits, the latter being rejected 
at > 99.99 CL. However, this is somewhat driven by our 
choice to correct the linestrengths to obtain a given age in 



the highest velocity dispersion bin, particularly for H/3. 
If instead, we attempt to fit H/J in all velocity dispersion 
bins simultaneously, then EXP and FR remain poorer 
fits than SSP but are no longer formally rejected. The 
OSP model, however, is still rejected at a high confidence 
level. 




log(CT) log((7) 



Fig. 3. — Linestrengths as a function of velocity dispersion. The 
Shapley binned linestrength data are shown by the black lines and 
circles, with error bars representing the uncertainty in the median. 
The predicted line index values for the best fit model parameters 
for each model are also shown: OSP (red short - long dashed); 
SSP (green dotted); EXP (blue short dashed line); AQ (cyan long 
dashed); FR (magenta dot - short dashed), STR orange dot - long 
dashed). Note that the OSP model is offset even in the highest 
velocity dispersion bin because we have adopted an age of 13 Gyr 
for this model rather than the default of 10.8 Gyr. It is clear that, 
with the exception of OSP, the models are nearly degenerate in 
most line indices. 



4.1.2. Fits to Scatter 

With at least three lines, we can iteratively calculate 
the best fit intrinsic scatter S'1.2,3 in the three model 
parameters A, 2, 3, where where Pi represents the age- 
related parameter (log(tssp) for the SSP case) and P2,3 
are metallicity and a-enhancement respectively: 




where the derivatives -^p- are calculated numerically from 
the model grids. 

The above formula assumes that, at a given a, the 
intrinsic scatter in age, for example, is independent of 
that of metallicity. However, we kno w that these pa - 
rameters are not in fact ind ependent. IWorthevI (|1994l ). 
and later iTrager et al.l (|2000D showed that there is a cor- 
relation between ag e and metallicity at a fixed mass. 
'Smith ct al. (2008a) found d[Z/ H]/dlog{tssp) = -0.68 
and 9[a/Fc]/91og(tssp) — 0.34, using the same Shapley 
data as in this paper. 
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TABLE 2 

Best-Fit Model Linestrength Indices by Bin and Model 



BkJ I Li^e Median (A) McdZj" ^ SSP AQ STR ixP FR OSP 



(T in km/s 


Observed Data 


Model Best Fits 


irtange of a 


HAT? 


U- 1 


U.UoO 


U.ZO 1 


U.DOU 


U.UOO 


U.DD ( 


n (^fi7 


u.u 1 y 


u. ozo 


= 35 ^ 59; 




-U.ooD 


U.UDO 


U.zuo 


-U. ( ( 


n 7fifi 

-U. (DO 


n Tfio 

-U. ( DZ 


-U. ( 00 


n 7/1 fi 

-U. ( 4o 


1 nfi/1 

-l.Uo4 






/I /I 1 v 


U. 144 


U.OIU 


A /I no 

4.4UZ 


A QOV 

4.oy / 


A /1 1 

4.4iy 


A A fin 

4.4DU 


A AA'V 
4.44 ( 


'i 09c; 
vS.yzo 


= 1.706 






U.UDU 


u.zzz 


i.y4D 


i.yoz 


1 01 7 

i.yi / 


1 fifl7 
l.oD ( 


l.oDD 


l.Oo ( 




r eoUlo 


O.U4i 


n 1 1 n 
U. IIU 


n /I Of; 
U.4UO 


O.Uoo 


o.UOo 


nc;n 
o.UoU 


c; no^; 

O.UZD 


nQ /I 
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A ■ 
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4.31 
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n /1 70 
U.4 / z 


= 59 ^ 84 


rl'yr 


-0.893 


0.102 


0.466 


-0.806 


-0.794 


-0.795 


-0.785 


-0.752 


-1.177 






4.Do4 


U. loo 


U.D ( 4 


A f^Ofi 

4.oyD 


A E^fiC 




y1 fi'i7 


A fifiQ 


4.UO0 
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1 01 
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1 Q01 
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l.D ( 1 




PGOUlO 


0. IZO 


n 1 1 1 

U. 11 1 


u.4oy 


99fi 


0. ZZO 


0. ZUO 


0. IDO 


^ 99*4 
0. ZZo 


4. 1 00 


Ngal = 37 


Mp-h^l 77 


3.761 


0.068 


0.300 


3.739 


3.739 


3.736 


3.728 


3.750 


3.738 










X ■ 


6.64 


7.91 


9.67 


16.5 


15.1 
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Range of cr 


IT ATT 


n /I on 

u.4zy 


U.U4Z 


U. 10v5 


n /I Qi 
U.4ol 


U.4o4 


U.4oD 


n /I /I /I 
U.444 


U.4oo 


n OC7 
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= 84 ^ 124; 


rl'yr 




U.U40 


n 1 QA 

U. lo4 


-1. lyo 


1 1 on 
-1. lyu 


1 1 on 
-1. lyu 


1 1 fi/1 
-1 . lo4 


1 1 flO 

-1. IDy 


1 /I 'i'i 


(log(<T)> 




o.UUo 


u.uyu 


U.oo 1 


4.yoD 


A Oflfi 
4.yD0 


A OVI 

4.y ( 1 


/I 001 

4.yy 1 


/I Ofifi 

4.yoo 


A A fiO 


= 2.022 


T T /D 

Hp 
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U.UvSo 


U. 104 


1 C01 


1 QClQ 

l.oUo 


1 70 7 
1. < y < 


1 'yt^ti 
1. (OD 


1 '7t^A 

1. /o4 


1 fioc; 
l.DZO 




peouiD 




U.Uoo 


U.oDZ 


/lOfi 

o.4yo 


/I fit; 
0.4oO 


/1 79 
0.4 / Z 


0.4o / 


/I /I fi 
0.44o 


r\A ^ 

0.U40 


Ng,l = 37 


Mp^h'Sl 77 


4.153 


0.070 


0.329 


4.132 


4.132 


4.132 


4.128 


4.127 


4.054 










v^- 

A ■ 


6.49 


7.55 


8.70 


15.0 


16.2 


136 


R ange of cr 


H5F 


0.354 


0.019 


0.166 


0.346 


0.347 


0.347 


0.348 


0.346 


0.281 


= 124 171; 


H7F 


-1.464 


0.046 


0.200 


-1.409 


-1.401 


-1.401 


-1.399 


-1.384 


-1.580 


{log(a)> 


Fe4383 


5.018 


0.082 


0.344 


4.972 


4.978 


4.980 


4.990 


4.992 


4.557 


= 2.160 


H/3 


1.792 


0.042 


0.192 


1.725 


1.723 


1.718 


1.688 


1.678 


1.597 




Fe5015 


5.446 


0.071 


0.307 


5.480 


5.501 


5.496 


5.473 


5.488 


5.216 


Ngal = 37 


Mgb5177 


4.363 


0.066 


0.309 


4.382 


4.377 


4.377 


4.375 


4.377 


4.533 












4.64 


5.51 


5.80 


8.28 


10.9 


89.9 


Range of cr 


HfiF 


0.198 


0.042 


0.181 


0.198 


0.199 


0.200 


0.201 


0.200 


0.124 


= 171 314; 


H7F 


-1.689 


0.052 


0.240 


-1.689 


-1.687 


-1.686 


-1.682 


-1.673 


-1.794 


(log(<x)> 


Fe4383 


5.219 


0.067 


0.290 


5.219 


5.222 


5.220 


5.227 


5.225 


5.116 


= 2.312 


H/3 


1.613 


0.039 


0.179 


1.613 


1.611 


1.609 


1.601 


1.584 


1.537 




Fe5015 


5.671 


0.068 


0.305 


5.671 


5.669 


5.667 


5.666 


5.668 


5.569 


Ng,i = 36 


Mgb5177 


4.692 


0.054 


0.250 


4.692 


4.693 


4.693 


4.692 


4.690 


4.714 










X': 


0.00 


0.00 


0.02 


0.14 


0.63 


15.8 








Bins 1 


4, Total x^: 


20.91 


24.8 


28.5 


47.2 


50.1 


409 



4.1.3. Distribution of Stellar Population Parameters as a 
Function of Velocity Dispersion 

Having outlined how the line index scatters are pre- 
dicted given correlated scatter in the model parameters, 
we proceed to fit these scatters to the observed values. 
Table |3] presents the model parameters yielding the best 
fits to the line index data. Also tabulated are the fit scat- 
ters in each parameter for each bin and model. Figures^] 
and [5] illustrate the parameter versus velocity dispersion 
relations, and their associated scatter. We also show the 
linear fits to the parameters, ignoring the lowest a bin, 
i.e. using data with a > 70 kms~^ (log(cr) > 1.844); the 
quoted zero points are calculated at log(cr) = 2.00. 

4.2. Systematics 
4.2.1. Aperture Effects 

Our data are observed through a 1" diameter fiber, cor- 
responding to 0.95 kpc radius at the distance of Shapley. 
Thus the stellar populations found above correspond to 
the stars within that aperture. Nevertheless, for some 
purposes, it is interesting to consider how this might be 
extrapolated to refer to all stars within the effective ra- 
dius. Re- As our data are all at a common distance we 
cannot solve for population gradients internally. Internal 
radial gradients in early-type galaxies have been studied 
by a number of authors, and are summarized in Table 4 
of iRawle et al.l (|2008D . They find that gradients in age 



To allow for correlated scatter, we modify equation 
([3]) as follows. We allow for a scatter in the age pa- 
rameter, and a correlated scatter in metallicity and a- 
enhancement. The sense of the correlation is that if a 
galaxy is older than the median by 0.1 dex, it is more 
metal-poor by -0.068 dex. We also allow for "extra" or 
residual uncorrelated scatter in the latter two parame- 
ters. We can now predict the scatter in each index as 
a function of the predicted scatter in each (SSP) model 
parameter as 




where t represents the age-related parameter (e.g. t = 
log(tssp) for the SSP case) and P2,3 are metallicity and 
a-enhancement respectively. 

By examining the response of the indices in the other 
star formation history models, in comparison with the 
SSP, it is possible to calculate the slope of the corre- 
lation between, for example, metallicity and the "age" 
parameter t (e.g. t in the case of EXP models). These 
correlations for both metallicity and a-enhancement for 
all models are tabulated in Table [H 
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TABLE 3 
Best-Fit Model Parameters by Bin 



Range of 


cr (km/s) 






(Age-) (Gyr) 




([a/Fe]) 






S([Z/H])res 


5([Q:/Fe])rcs 


Single Stellar Population Model (SSP); Age = time since single burst 


35 - 


59 


1.706 


3.14 


5.91 


-0.002 


0.189 


25.6 


0.186 


0.001 


0.105 


59 — 


84 


1.844 


6.64 


5.53 


0.089 


0.186 


7.32 


0.212 


0.119 


0.001 


— 




9 099 


6.49 


7.08 


0.185 


0.214 


19.7 


0.142 


u.u / y 


U.U^D 




1 71 


9 1 fin 


4.64 


9.10 


0.190 


0.260 


21.6 


0.151 


u.u / 


n OAS 


171 — 


314 


2.312 


0.00 


10.8 


0.240 


0.280 


33.4 


0.141 


0.070 


0.041 








Abruptly Quenched Model (AQ); A 


ge = time since quenching 






35 — 


59 


1.706 


3.84 


3.17 


-0.009 


0.183 


39.9 


0.228 


0.001 


0.115 


59 - 


84 


1.844 


7.91 


2.66 


0.081 


0.180 


8.61 


0.331 


0.121 


0.001 


84 — 


124 


2.022 


7.55 


3.99 


0.180 


0.211 


22.1 


0.238 


0.085 


0.039 


124 - 


171 


2.160 


5.51 


6.15 


0.188 


0.257 


24.3 


0.271 


0.080 


0.047 


171 - 


314 


2.312 


0.00 


9.06 


0.240 


0.280 


44.9 


0.136 


0.090 


0.037 








Strangulation Model 


(STR); Age 


= time since quenching 






35 - 


59 


1.706 


4.31 


5.44 


-0.018 


0.175 


41.2 


0.123 


0.001 


0.116 


59 - 


84 


1.844 


9.67 


5.05 


0.056 


0.169 


4.95 


0.218 


0.109 


0.001 


84 - 


124 


2.022 


8.70 


6.14 


0.169 


0.206 


25.5 


0.087 


0.084 


0.001 


124 - 


171 


2.160 


5.80 


7.90 


0.184 


0.255 


22.2 


0.127 


0.074 


0.001 


171 - 


314 


2.312 


0.02 


10.57 


0.238 


0.280 


43.6 


0.117 


0.063 


0.036 








Exponential Star I'ormation Rate Mode 


(EXP); 


i-ge = e- 


folding time 






35 - 


59 


1.706 


7.44 


2.28 


-0.040 


0.161 


27.2 


0.711 


0.001 


0.095 


59 - 


84 


1.844 


16.5 


2.39 


0.005 


0.145 


18.6 


0.829 


0.109 


0.001 


84 - 


124 


2.022 


15.0 


2.10 


0.133 


0.190 


30.0 


0.366 


0.098 


0.039 


124 - 


171 


2.160 


8.28 


1.81 


0.160 


0.246 


27.5 


0.613 


0.085 


0.015 


171 - 


314 


2.312 


0.14 


1.44 


0.233 


0.277 


15.2 


0.903 


0.082 


0.001 


Frosting Model (FR); Age = time since secondary burst 


35 - 


59 


1.706 


7.88 


1.37 


-0.055 


0.160 


32.8 


0.146 


0.001 


0.105 


59 - 


84 


1.844 


15.1 


1.23 


0.000 


0.132 


13.9 


0.159 


0.122 


0.001 


84 - 


124 


2.022 


16.2 


1.48 


0.119 


0.184 


27.1 


0.085 


0.099 


0.001 


124 - 


171 


2.160 


10.9 


1.78 


0.153 


0.241 


26.0 


0.097 


0.099 


0.001 


171 - 


314 


2.312 


0.63 


2.49 


0.224 


0.274 


19.9 


0.199 


0.087 


0.001 


Old S 


ngle Stellar Population Model (OSP); A 


ge = time since single burst (13Gyr) 




35 - 


59 


1.706 


75.6 


13.00 


-0.194 


0.303 


75.4 


0.000 


0.050 


0.133 


59 - 


84 


1.844 


107 


13.00 


-0.151 


0.303 


47.0 


0.000 


0.137 


0.067 


84- 


124 


2.022 


136 


13.00 


-0.059 


0.272 


35.6 


0.000 


0.109 


0.05 


124 - 


171 


2.160 


89.9 


13.00 


0.044 


0.338 


38.4 


0.000 


0.105 


0.057 


171 - 


314 


2.312 


15.8 


13.00 


0.173 


0.293 


50.2 


0.000 


0.096 


0.041 



EXP: t; AQ: (aq; FR: tpa; STR: isTR*" Age Parameters: SSP: log(tssp); EXP: r; AQ: log(tAQ); FR: 



Age Parameters: SSP: tssp; 
log(tpR); STR: log(tsTR) 



TABLE 4 
Parameter Correlations 



Model 


a[z/H] 
ai° 


S[a/Fe] 


SSP 


-0.68=^ 


0.34=^ 


EXP 


0.172 


-0.086 


AQ 


-0.367 


0.183 


FR 


-0.523 


0.262 


STR 


-0.610 


0.305 



ISmith et al.l ll2008al') '^ t = 
model age parameter. See ta- 
ble |3] note b. 

and a-enhancement are small, but the gradient in Z is 



AZ/Alogra 



-0.18 ± 0.02. We can correct our data 



to i?e using this result, together with eq. [H We find 

\ZlH\{R:)-\ZlH\{av)^ 

-0.043 log(cr) < 2.23, 

-0.50(logCT/2.23) - 0.043 \og(a) > 2.23. 



(5) 



Thus the correction is small: —0.08 for the highest ve- 
locity dispersion bin, and ^ —0.04 for the other bins. 

4.2.2. Line Indices: Choice of Index and Effect of Index 
Calibration Corrections 

It is interesting to see whether the conclusions, par- 
ticularly those for age, change significantly depending on 
which lines are used in the analysis. Here we focus on the 



Balmer lines, reproducing the fits with only one Balmer 
line rather than all three simultaneously. The results of 
these tests are shown in Fig. [6l We find that the median 
age in the second- lowest bin (logcr = 1.844) changes by 
less than 1 Gyr (15%) depending on which Balmer line is 
used. In H7 and H/3 we see a flattening or an upturn in 
age at the lowest velocity dispersion bin (as in Fig. [6]), 
but if we only use H(5, the lowest bin is younger. 

One concern regarding the Balmer line indices is the 
possibi lity of weak infillin g by AGN emission at high 
mass. ISmith et alJ (I200I Fi gure 7 shows the fraction 
of RSGs with AGN-like emission in the Shapley sample. 
At Mji < —20.5, and thus high mass and velocity dis- 
persion, there is still significant AGN contribution on the 
red-sequence. While galaxies with clear emission have 
been removed from our sample, it is possible that weak 
AGN emission may cause infilling in Hf3, particularly for 
the highest velocity dispersion galaxies. If so, it is possi- 
ble we have overcorrected HP by the procedure adopted 
in Section m An alternative calibration would have been 
to fit a zero-point iteratively to the average of all bins 
rather than to the highest velocity dispersion bin alone. 
In fact, if we fit the observed H(3 indices to all veloc- 
ity dispersion bins, then all models agree better with the 
predicted values of Hf3. The weak effect of this alter- 
native correction on the derived ages is shown in Figure 

El 
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slope = 0.662, z p. -0.842, r.m s -0 167 
I I I I I I I I I I I I I 

1.8 2 2.2 2.4 

log(a) 



2 2.2 2.4 
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log(CT) 



2 2.2 
log((7) 



2 2.2 
log(CT) 
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log((7) 
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log((7) 



2 2.2 
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Fig. 4. — Best-fit model parameters as a function of \og{a). The filled points show the median values of the parameters and the solid 
black lines show the itlo" range. Fitted slopes (shown by red dotted lines) are to data with log((T) > 1.844 only. Rows from top to bottom 
correspond to SSP, AQ, and STR. Age parameters are as in Table [3] Note increasing metallicity, o-enhancement and age with increasing 
velocity dispersion. Synthetic cluster density contours are shown in grey, at O.Tdex increments (see Section[5]l. 

A related concern is the effect o n our conclu s ions o f 
the choice to calibrate all lines to iNelan et al.l (|2005f ). 
which itself was not calibrated to the Lick system. To 
test the effect of alternative calibration schemes, we 
have also compared our Shapley linestrength-cr relations 
to those from 12 Com a early-type galaxies observed 
by iTrager et al.l ()2008D . which we expect to be well- 
calibrated to the Lick system. We then fit additional 
offsets to our Nelan-calibrated Shapley linestrength-cr re- 



lations to match the lTrager et al.l ()2008D linestrengths as 
a function of a. In particular, the Balmer line offsets were 
estimated as HJF H7F,H/3 ~ ( ).12,0 2,0.1 (in addition to 
the correction to INelan et al.l (|2005f ) tabulated in Table 
[T|). We then refit our models with this alternative cali- 
bration scheme, and find that it yields ages that are only 
slightly younger (by < 0.5 Gyr or 10% for the youngest 
galaxies). This change in stellar population parameters 
is smaller than other systematics, such as the effect of 
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jlope-0.000, z. p. -1.114, r.m.s. -0.000 
I I I I I I I I I I I I I 
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Fig. 5. — Best-fit model parameter trends with log{(7). Rows from top to bottom correspond to EXP, FR, and OSP. Lines, contours and 
symbols are as in Figure |4] Age parameters are as in Table |3] Note increasing metallicity, a-enliancement and age with increasing velocity 
dispersion. 

using only Hf3 as discussed in the previous paragraph. 



4.2.3. Systematics Associated with Model Parameters 

In Section [3l in certain models, star formation was 
assumed to start 13 Gyr ago, which corresponds to a 
rather high redshift. As an alternative, we will consider 
a different starting time and investigate its effect on the 
fits. For this purpose, we choose a redshift of 2, corre- 
sponding to approximately 10.3 Gyr, as a more realistic 
choice, as this is close to the peak star formation in the 



universe ms fc Beacoml[2006l) . When we refit our 

star formation history models, we find that the CSFH 
age parameters shift towards older ages to compensate 
for the absence of very old stars. For example, for the 
EXP model, for the {a) = 70 kms~^ bin, r is reduced, 
from its default value of 2.4 Gyr, to only 1.7 Gyr. In 
the AQ model, the lookback time to quenching increased 
from 2.7 Gyr to 3.3 Gyr. Other velocity dispersion bins 
are affected similarly. We also find that metallicity in- 
creases by approximately 7% in the "mostly old" models. 
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H6F only 
H7F only 

only 
HdF, H7F, 



1.6 



1.8 



2.2 



2.4 



log(CT) 



Fig. 6. — SSP age as a function of log{a) for different clioices of 
Balmer line used in the fits. (The Fe and Mg lines are used in all 
cases.) In the H/3* case, we have fit H/3 usi ng a calibration fit to all 
velocity dispersion bins (see Section l4.2.2l for details). In all cases, 
the choice of Balmer index makes little dilference at the youngest 
ages. 

whereas the intermediate age models show no significant 
change in metallicity. Finally a-enhancements remain 
virtually unchanged from the 13 Gyr case. We will dis- 
cuss the effects of these changes on stellar mass-to-light 
ratios in Section [7TT1 

4.2.4. Effect of Other Population Synthesis Models 

iSmith et al.1 (|2009b[ ) found a generally good agreement 
between o ur TMBK - based SSP models and the EZ-Ages 
model of ISchiavonI (|2007D for age, metallicity and a- 
enhancement. However, they noted that there was a 
small [Fe/H]-dependent bias in the derived ages in the 
sense that the EZ-Ages ages were younger by ~ 0.25 dex 
per dex of [Fe/H]. Given the small range in ([Fe/H]) cov- 
ered by our sample (only 0.05 dex for a > 70km s~ ), the 
effect of the choice of models on the Age-cr relation and 
hence on colors and mass-to-light ratios is negligible. 

4.3. Comparison with previous work 

For all models, the general stellar population trends 
observed here for a > 70 kms~^, are that age, metal- 
licity and a-enhanccmcnt all increase with velocity dis- 
persion, are in agreement with a number of previ- 
ous results ([Caldwell et al.l l200l iHeayens et al.l 120041: 
Nelan et al.l 120051: iThomas et al.1 120051: iGallazzi et all 
2006; Smi th et alhoOTl) . 

The scaling relations we find for a > 
tabulated in Table [51 Note that these 



70 



kms~ are 
are based on 

the fits to synthetic clusters (discussed below in Sec- 
tion [5]) and are weighte d to low velocity d ispersions. 
We also show results from lSmith et al.l ()2007l ). using the 
same Shapley data (and indices) for all galaxies and for 
cr > 100 kms^^. In general our scaling relations are 
between the "all a" and "(t > 100 kms~^" relations 
found bv ISmith et al.l ()2007[) . although our scaling re- 
lation for [Z/H] is somewhat flatter. In comparison with 



I Nelan et al.l (|2005f ). the scaling of age with velocity dis- 
persion is similar, however this analysis yields shallower 
relationships with [Z/H] and [a/Fe]. Note that our de- 
termination of a-enhancement is based on Mgb and so 



Mg/Fe] abundanc e ratio, whereas 
Nelan et al.l (|200l was also based 



primarily reflects the 

the a-enhancement of ^ ^ 

on CN . A full table of comparisons is given in lSmith et al.l 
(|200l . 

The lowest-velocity dispersion bin (cr < 70km s^^) 
shows a either a flattening or an increase in the mean 
age of galaxies, in comparison to cr > 70km s~^. This is 
as a result of the downturn or flattening of the Balmer- 
line indices seen in Fig[31 We note that this change in the 
ages o f low velocity dispe rsion galaxies can be seen in the 
fits of Smi th et al.l (|2007f ). in the sense that their scaling 
relations are flatter when all galaxies are included and be- 
come steeper when only high velocity dispersion galaxies 
(cr > 100) kms~^ are used in the fit. Furthermore, this 
trend is also seen in central Coma cluste r galaxies at low 
velocity dispersions (S mith et alll2009bD . Consequently, 
we argue that this effect is real, and signals the end of 
"downsizing" on the cluster red-sequence. We discuss 
selection effects that may affect the derived stellar popu- 
lation parameters in Section but conclude there that 
correction for such effects would strengthen our conclu- 
sion that this "upsizing" effect is real. 

In Section we also show that the "pinching" of the 
scatter S in the lowest velocity dispersion bin is a selec- 
tion effect. In reality, the scatter in age increases at de- 
creasing velocity dispersion. This increased scatter may 
also help to reconcile previous ly contradictory results 
in the literature. For example, ISanchez-Blazguez et al.l 
(|2006f ) find no "downsizing" effect in the Coma cluster 
but do find downsizing in Virgo. In both clusters they 
have a number of objects with high-cr plus a small num- 
ber (4-5) with log(cr) < 1.844. Given the large intrin- 
sic scatter in age that we find at low-cr (>0.2 dex or 
~ 3 Gyr), we argue that a sample of 4-5 galaxies is too 
small to secure a robust detection of the me an age of the 
low-cr population. Similarly, with regard to lTrager et al.l 
(j2008), who also find no downsizing with a 12 galaxy 
sample, we argue that a sample of this size is too small 
to detect the downsizing effect, given the large intrinsic 
scatter in the age-a relation. 

5. SYNTHETIC GALAXY CLUSTERS 

In section|4l we have shown how we derive the medians 
and dispersion of age, metallicity and a-enhancement as 
a function of velocity dispersion. For a given choice of 
stellar population parameters, P, we can generate stellar 
mass-to-light ratios, and hence colors and mass-to-light 
ratios, as outlined below in Section 15.21 To populate a 
synthetic cluster, the additional ingredients required are 
the velocity dispersion distribution, effective radii and 
stellar masses. 

5.1. The Velocity Dispersion Distribution 

The next step is to generate a synthetic cluster pop- 
ulated by galaxies with appropriate statistical distribu- 
tions of velocity dispersion, effective radius, stellar mass, 
age and metallicity, as well as derived parameters (lumi- 
nosity). Because velocity dispersion is the driving param- 
eter, a prerequisite input is the velocity dispersion distri- 
bution (VDD). We adopt a lower limit of cr = 20km s~^ 
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TABLE 5 
SSP Scaling Relations 



Reference 1 




a range 


Age 


Z/H 


a/Fe 


Nelan et al. 


2005')! 


50km s 


1 < cr < 400km 


0.59 ±0.13 


0.53 ± 0.08 


0.31 ±0.06 


Smith et al. 


^2007) 


30km s~ 


1 < o- < 300km 


0.52 ±0.06 


0.34 ± 0.04 


0.23 ± 0.04 


Smith et al. 


f2007) 


100km s■ 


< o- < 300km s-l 


0.64 ±0.12 


0.38 ±0.09 


0.36 ± 0.07 


This work 




70kms- 


1 < (T < 300km 


0.66 


0.27 


0.22 



for the velocity dispersion distribution for our synthetic 
clusters. While we will not create velocity-dispersion- 
limited samples with velocity dispersions this low, be- 
cause of the large scatter in the Faber- Jackson relation, 
it is necessary to include low velocity-dispersion galaxies 
in order to ensure completeness at a given magnitude. 

To mode l the VDD, we begin with the fits of 
ISheth et al.l ()2003[ ) to a sample of de Vaucouleurs-profile 
early-type galaxies in the SDSS 



(a)da — 



exp 


- (a/a,)^" 




r (a//?) 



/3 



da 



(6) 



where the best fit parameters are 0* = 0.0020 ± 0.0001, 
a = 6.5±1.0, 15 = 1.93±0.22ando -., = 88.8±17.7kms-i. 

According to lSheth et al.l ()2003l . Figure 6), one expects 
that late- type galaxies dominate the VDD, at least for 
velocity dispersions lower than 200 km/s. However, the 
exact level of the contribution of these galaxies to the 
VDD is very uncertain. Nevertheless, to the velocity dis- 
persion limit to which we would like to synthesize galax- 
ies (~ 20km s~^), we expect non-early- types to make a 
non-negligible contribution to the RSG cluster popula- 
tion. Thus we model the "faint" end of the VDD as a 
power law in a: 

0(cr) oc cr^ (o- < 200km s-i) . (7) 

We choose to fit the "faint-end" slope, ^, by adjusting 
it to fit the observed z ~ red-sequenc e dwarf-to-giant 
galaxy ratio (DGR) ()Gilbank fc Baloghl [20081. We find 
this slope to be ^ = —1.275, yielding a DGR of 3-4 at 
z = 0. The final adopted distribution of a is shown in 
Figure!?] For each of our CSFH models, we create 50000 
galaxies from this distribution, with cr > 20 kms~^. 

Thus the procedure is as follows: generate a velocity 
dispersion from the VDD described above. For a given 
velocity dispersion, the mean and scatter in the stellar 
population parameters are known (the binned values are 
interpolated or extrapolated where necessary). For each 
of the stellar population parameters, we can thus gen- 
erate a realization consistent with the statistical distri- 
butions, allowing for correlation between offset in "age" 
and the offset in metallicity, and likewise between "age" 
and a-enhancement. The distribution of the simulated 
populations are shown by the contours in Figures [4] and 

El 

5.2. Mass-To-Light Ratios and Colors 

Having generated the stellar population parameters for 
our synthetic, clusters, we now turn to predicti ng observ- 
ablcs, such as color. The models of Maraston (jMarastonl 
[1998, 2005) provide grids of stellar mass-to-hght ratios 
in both the Johnson-Cousins and SDSS filter sets, which 
can be used to compute colors as a function of age and 
metalUcity. Here we adopt the lKroupal (|2001h IMF with 



an 



-3 



1.5 



2.5 



log(c7) 



Fig. 7. — The normalized velocity dispersion distribution (VDD) 
of early- types from Shcth ct al. (2003), with extrapolation to low 
velocity dispersions shown by the dashed line. Late-type galaxies 
dominate at cr < 200km s~^. We approximate this with a power 
law (solid line) and fit the slope such that the resulting dwarf-to- 
giant ratio of all models reasonably matches the z observed 
dwarf-to-giant ratio of Gilbank & Balogh (2008). 

red/intermediate horizontal branch. Using the Salpeter 
IMF would increase the M^/L ratios by a factor of ~ 1.6 
([Thomas et al.ll2005D : we discuss the effects of the choice 
of IMF in Section [II 

We then fit scaling relations of mass-to-light ratio and 
color as a function of velocity dispersion. For the fits, 
only galaxies with velocity dispersion log(cr) > 1.844 are 
used, as these are not affected by the "upturn" in age at 
low velocity dispersion (Figures [H [Sj) . In the Appendix, 
Tables [H and [21 summarize the slopes and intercepts as a 
function of velocity dispersion for Johnson-Cousins and 
SDSS passbands. 

5.3. Effective Radii, Stellar Masses and Luminosities 

In order to predict luminosities and the Fundamental 
Plane, we need to assign a stellar mass and an effective 
radius to each synthetic galaxy. 

We assign an Re based on the a, following the broken- 
power-law model of our Coma data in Section [2. 21 allow- 
ing for scatter in this relation by generating a random 
Gaussian deviate with the same root-mean-square scat- 
ter as given in Figure [H Total mass is then assigned 
according to the virial relation 



M = C2 



cr'tRe 
G 



(8) 
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where here we assume Ce ~ cr (see Section |6^ . The fac- 
tor C2 depends on the himinosity distribution, velocity 
dispersion anisotropy a nd rotation. Perhaps surprisingly, 
ICappellari et al.l (|2006f) found that a constant C2 = 5 was 
an excellent approximation to detailed dynamical mass 
reconstructions. This indicates that there is little evi- 
dence for non-homology, i.e. a dependence of C2 on mass 
or velocity dispersion. The same result, that early- types 
are homologous, has also bee n found using gravi tational 
lensi ng as a mass e stimator ()Bolton et al.l [20081) . Note 
that ICappellari et~al, ( 200Q) measured Re in the /-band, 
which is close to the r band used here for measuring i?e, 
so their value of C2 may be transferred to this analysis. 
Denoting the fraction of the total mass in stars as /*, 
we then have M, = /*M. Clearly for our analysis, C2 
and /, are degenerate. We adopt C2 = 5 and /* = 1 (no 
dark matter), and return to discuss these assumptions in 
Section O 

Finally, the magnitudes are then given by, for example 
in the i?-band, 

M. = Af.,-2.51og(^) (9) 

where the stellar M^jLji is given by the Maraston models 
for the known stellar population parameters, constructed 
as described in Section [31 

5.4. Selection Effects at Faint Magnitudes 

As noted above, the model fits show an increase 
in all ag e-related param eters at low velocity disper- 
sion. Kels on et al.l (|2006[ ) pointed out that a magnitude- 
limited sample will lead to biases in stellar population pa- 
rameters at fixed mass because, for example, older galax- 
ies will be fainter and hence will be exclude d by the mag- 
nitude limit. This bias was also studied bv lGraves et al.l 
{2007), who calibrated the bias empirically, and found a 
larger correction to metallicity than to age. 

Here, we use our synthetic cluster populations to de- 
termine the effect of this bias on our modeled stellar 
population parameters. Our Shapley sample was lim- 
ited to TUfi < 18, or Mfl < —18.65, but our synthetic 
clusters are limited by velocity dispersion, not by mag- 
nitude. By removing the synthetic galaxies fainter than 
— 18.65, we can quantify the biases that may result from 
this selection by refitting the model parameters to the 
magnitude-limited sample. Figure |5| illustrates the ef- 
fect, in the case of the SSP model. We find that, when 
faint galaxies are excluded, the remaining low-cr galax- 
ies are younger, more metal-rich and less-a enhanced. 
One might naively expect the bias in metallicity to go in 
the opposite sense: i.e. metal-poor galaxies are brighter 
and hence they should be preferentially included in a 
magnitude-limited sample. However, in our SSP model, 
the scatter in age and metallicity at fixed a is assumed 
to be anti-correlated (Section I4.1.2p . Since the effect of 
age on mass-to-light is greater than that of metallicity, 
the bias towards younger brighter galaxies implies that 
we are also selecting relatively metal-rich galaxies. The 
correlated scatter between age and a-enhancement goes 
in the opposite sense, and hence so does the correction 
to a-enhancement. The effects of this selection bias are 
small, however, and affect only the lowest velocity dis- 
persion bin. Now because the original Shapley data were 
indeed magnitude-limited, then to account properly for 



this selection effect, a correction should be applied (in 
the opposite sense to that found when faint galaxies are 
removed from the synthetic cluster sample). The cor- 
rected stellar population parameters are shown by the 
dashed line in Figure HI Thus the correction, if applied, 
would make the low velocity dispersion galaxies older, 
less metal-rich and more a-enhanced. The same bias also 
affects the scatter in stellar population parameters: be- 
cause galaxies are removed from the magnitude-limited 
sample, the scatter in, for example, age is reduced. Thus 
a correction for selection would increase the scatter at low 
velocity dispersion. We do not correct for these effects 
in this paper, because Figure |8| shows that the effects are 
small. Nevertheless, the sense of the correction strength- 
ens our conclusion that the "upsizing" at low velocity 
dispersion is real. Note, however, that the selection does 
not affect the robust "downsizing" trend seen at higher 
velocity dispersion (cr > 70kms~^). 

6. RESULTS 

For each model, we now have stellar population pa- 
rameters and derived parameters, such as magnitudes 
and colors, for a synthetic cluster galaxy population. We 
can now compare these to independent observations in 
order to constrain the model. Later in this section, we 
will compare the predictions to scaling relations such as 
the FP. We begin however with a comparison sensitive 
only to the stellar populations: the color-cr relation. 

6.1. Comparison of color-a relations 

Figure |9| compares predicted u — g,g — i and i--K colors 
in our models to the Coma cluster data measured in an 
aperture of the same physical size (~ 0.95 kpc, 2") as the 
Shapley line index data. The quantitative comparisons 
are restricted to galaxies with cluster-centric radii ^200/2 
to match our models derived from the central regions of 
Shapley clusters. 

Although variable a-enhancem ent is used to fit the 
line indices using the models of iThomas et al. (l2003l . 
I2004D . the stellar M*/L predictions of lMarastonI ()2005D 
do not allow for the effect of a-enhancement. We can 
estim ate the effect s of g -enhancement using the mod- 
els of ICoelho et al.l (|2007D . who compute tracks and col- 
ors for an a-enhanced model as well as models with so- 
lar abundance ratios. Their colors are tabulated as a 
function of [Fe/H] and [a/Fe]. Fitting the colors as a 
function of age, a-enhancement, and metallicity (where 
[Z/H] = [Fe/H] + 0.75[a/Fe] f or the enhancement pat- 
tern adopted bv ICoelho et al.l (|2007l )). we can find the 
differential effect of [a/Fe] (at fixed [Z/H]) on the colors. 
Since a-enhancement is a weak function of velocity dis- 
persion, in principle both the slopes and the zero-points 
of the color-cr relations could change. Using their self- 
consistent SSP model in a differential sense, ie. to derive 
an a-enhancement-dependent color correction, we find 
that the optical colors shift blueward with increasing a- 
enhancement, at fixed Z/H. We illustrate this for the SSP 
case for u — g and g — i 'u\ Figure |9l where the lower SSP 
model is with a-enhancement. 

Note that while our models do not allow for dust, in 
reality, of course dust may be present, particularly for 
low-cr galaxies, that may be of later morphological type. 
Correcting the data for this possible effect would shift 
the low-cr data blueward, and would make slope steeper. 
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Fig. 8. — Effect of a magnitude limit on the modeled SSP stellar population parameters. The Shapley data, on which the stellar 
populations of our synthetic clusters are based, were limited to M/j < —18.65. Thus the derived stellar populations, and hence the mock 
clusters, may be affected by biases due to this magnitude limit. However, having fixed the stellar populations, there is no explicit magnitude 
limit in the synthetic clusters, i.e. some galaxies have magnitudes fainter than —18.65. This allows us to use the synthetic clusters to 
test the biases due to the magnitude limit. Black contours represent galaxies with M^j < —18.65, whereas red contours represent galaxies 
with Mji > —18.65 which would have been excluded had our synthetic clusters been magnitude limited. The dot - dashed red lines are 
fits to the stellar populations of the synthetic clusters with the faint galaxies included in the entire population, whereas the solid black 
lines are the fits when the faint galaxies are removed. When faint galaxies are removed from the sample, wo find that, at the low velocity 
dispersions, the population is generally younger, more metal and less o-enhanced. The scatter in each stellar population parameter also 
decreases. Therefore, in order to correct for this selection effect in the real data, we would have to reverse the offset from this test, resulting 
in an underlying population that is older, more metal-poor and more a-enhanced with greater scatter in all parameters at dispersions 
log{(T) < 1.85. The dashed black lines show our best estimate of the underlying population parameters after correction for this effect. 



In u — (7, without allowing for a- enhancement, all mod- 
els are too red. The best fit slope is the FR model and 
all the other models overpredict the observed slope, by 
factors as large as 33% (SSP). However, if we allow for a- 
enhancement, the SSP prediction for u—g shifts blueward 
by -0.090 at logcr = 1.84 and -0.136 at log cr = 2.31. 
This makes the SSP model very slightly bluer (by a few 
hundredths of a magnitude) than the Coma data, and 
flattens the slope to 0.45, in good agreement with the 
observed slope (0.43 ±0.04). We have not calculated the 
effect of a-enhancement on the other models, but one 
might expect the sense of this correction to be the same 
as for the SSP case, i.e. to flatten the model slope and 
shift the model to bluer colors on average. Note that u—g 
is, however, sensitive to dust, and so the small reddening 
of the data may be due to dust extinction. 

In the g — i color, all of the models are again too 
red and, with the exception of the SSP model, are also 
too shallow compared to the data. The effect of a- 
enhancement is weaker \m g — i than in u — g: the 
models shift blueward by —0.026 at loga = 1.84 and 
-0.039 at logCT = 2.31 for the SSP case. Thus allow- 
ing for a-enhancement, all models remain slightly too 
red (by ^ 0.04 mag in the SSP case), and of course, 
a disagreement in this direction cannot arise from dust. 
The SSP {g - i) - \oga slope flattens slightly to 0.25, 
but is still a reasonable fit to th e slope of the data 
(0.31 ±0.03). iMaraston et~aI1 ([2008fl have no ted that, by 
using empirical libraries from iPicklesI (jl998l ). the colors 
of solar-metallicity models change, particularly in g and 
r. However, the sense of the change would be to make 
the models even redder in g-i. In contrast to the good 
agreement with the a-enhancement-corrected SSP case, 
the "mostly-old" models (EXP, FR, OSP) have slopes 
only two-thirds of that observed before considering a- 
enhancement, which one expects would flatten them fur- 
ther. Furthermore, correcting for dust, if preferentially 



present in low-<7 galaxies, would steepen the observed 
slope, exacerbating the conflict with the "mostly-old" 
models. 

In i — if, again the model slopes are all too shallow. 
The OSP model has the strongest metallicity trend, and 
since i — K \s more metallicity sensitive than age sensi- 
tive, it yields the steepest i — K slope. However, the 
SSP model is the next ste epest model a fter t he OSP 
case. While the models of ICoelho et al.l (|2007[ ) do not 
pr edict the effec t of a -enhancement on j — if, those 
of iPercival et al.l ()2009l ) do predict its effect on I — K: 
~ 0.04 for the a-enhancement relevant at cr > 70 kms^^. 
Thus qualitatively we would expect a-enhancement to 
shift colors to the blue and very slightly flatten the pre- 
dicted color-cr relation, as with the optical colors. 

We conclude that, overall, the SSP model and, to a 
lesser extent, the AQ model are good fits to the opti- 
cal color-cr slopes, once the effects of a-enhancement are 
accounted for. The FR and EXP models do not fit the 
observed trends well in g — i. None of the models fit the 
observed i — K slope well, even the OSP model which 
has a strong dependence of metallicity on a. The com- 
parison of color-cr zero-points is less certain, particularly 
because of the effects of a-enhancement and of dust. No 
model fits the zero-point of, for example, g — i, although 
the SSP model zero-point is closest to the Coma data. 
We conclude that colors are not a strong discriminant of 
the CSFH models. 

6.2. The Faber- Jackson relation 

As described in Section 15.31 stellar masses are gen- 
erated based on the observed Re — <J relation, and from 
these and the mass-to-light ratios one can calculate mag- 
nitudes and surface brightnesses. 

Figure [10] illustrates the resulting Faber- Jackson (FJ) 
relations both for the Coma data and the models. 
Note that over the range of velocity dispersions cov- 
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Fig. 9. — Color - cr relations compared to model predictions. Colors are measured in a 2" aperture on Coma to match the 1" aperture of 
the Shapley data, at twice the distance. Fits to the Coma data are shown by the solid line, and robust fits to the medians, hence insensitive 
to outliers. Only galaxies within r20o/2 (filled circles) and cr > 70km s~^ are used in the fit. Line styles and colors as in Fig. \3\ except for 
the heavy green dotted line is the SSP model with a correction for Q-enhancement. The a-enhancement-correction reduces the color offsets 
in u — g and g — i, as discussed in the text. 



the observed slope of the FJ is close to 

^2 



5, 



ered here 

i.e. £ PC (7^. This is in agreement with previous work 
(|Matkovic fc Guzman|[2005i) . Only at large velocity dis- 
persions (log((T) > 2.23) does FJ approach the "classical" 
relation L oc a"^. This change in slope is a direct result 
of the change in slope of the log(i?e) — log(f7) relation 
discussed above (Section 15. 3p . The FJ relation is not a 
strong model discriminator because its slope and scat- 
ter are driven by the log(i?e) — logic) relation, and not 
by stellar populations and their associated mass-to-light 
ratios. 
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Fig. 10. — Coma Faber- Jackson relation compared to model 
predictions. Line styles and colors are as in Fig. |3] 
The fit to Coma data uses galaxies within a cluster-centric r2oo /2 
only and with cr > 70 kms~^. Note the transition in the data 
from L oc cr-^ to L oc cr'* at log cr ~ 2.2, resulting from break in the 
Re — cr relation (Figure [J}. The same transition is seen in FJ of 
the synthetic clusters, but here we simply plot the best linear fit 
to their FJ relation. 



6.3. The Fundamental Plane 

The FP is considerably tighter than the FJ relation, 
and hence offers a more powerful test of dynamical mass- 
to-light ratios. Figure [TT] shows the r-band FP data for 
Coma where we have chosen to regress the combination 

[l0g(i?e) - 0.32(^e)] on logcr. 
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Fig. 11. — The Fundamental Plane. The r-band data from Coma 
are reproduced in each panel by the open circles. The scatter in 
[log(i?e) — 0.32(/ie)] for Coma is 0.085. Line styles and colors are 
as in Fig. \3\ In order of increasing extra "tilt" relative to the data, 
these are SSP, AQ, STR, EXP, FR and OSP. We note however that 
the SSP slightly overproduces the scatter in the FP. However all 
models, with the exception of SSP and AQ predict, at fixed Re, 
too faint a surface brightness at low cr. 



If the tilt of the FP with respect to the virial relation 
were due only to variation in average stellar populations 
as a function of a, then we would expect the slope of the 
model FP to agree with that of the Coma data. Thus 
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the observed slope can be used to determine the addi- 
tional tilt that is due to DM and/or non-homology. We 
find that the SSP and AQ models require the least ad- 
ditional tilt from DM & non-homology. We discuss this 
in more detail in Section 17.31 below. Furthermore, the 
"mostly-old" (EXP, FR, OSP) models predict a surface 
brightness too faint compared to the data, at low ve- 
locity dispersions and at a fixed effective radius. The 
addition of dark matter would exacerbate this situation 
as it would increase the predicted velocity dispersion and 
so would move the models further to the right. 

An alternative way to address this issue is by con- 
structing the dynamical mass-to -light ratio directly 
from the observed FP parameters (iDressler et al.l 119871 : 
iBender et all Il992t iJorgensen et al l I1996D . These can 
then be compared to stellar mass-to-light ratio predic- 
tions, as shown in Fig. [T^l The most striking aspect 
of this is how closely the Coma dynamical M/L, as a 
function of cr, follows the distribution of stellar ages as 
a function of a (see Fig. |4]). Specifically we see the same 
steep slope at high a and the same upturn or fiattening, 
and increased scatter at low a. This strongly suggests 
that the FP tilt and scatter are driven primarily by stel- 
lar age eff'ects. 
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AQ slope = 0.455, z.p, = 0.413, r.m.s. = 0.063 

STR slope = 0.391, z.p. = 0.434, r.m.s. = 0.043 

EXP slope = 0.335, z.p. = 0.505, r.m.s. = 0.038 

FR slope = 0.174, z.p. = 0,561, r.m.s. = 0,031 

OSP slope = 0.196, z.p. = 0.563, r.m.s. = 0.032 
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Fig. 12. — The r-band mass-to-light ratio data from Coma. 
Observed data from the Coma Cluster are shown as black points, 
and the fit to these is indicated by the solid black line. Binned 
medians are plotted as grey squares. These bins represent the same 
log{(7) ranges used when fitting to Shapley (see Table |2j. Model 
fits to the stellar Mt/L are shown by the colored broken lines (line 
styles and colors as in Fig. I^. Again, we only fit these lines to 
log((T) > 1.844, see Section 15.41 We find the models consisting 
of intermediate mass stars (SSP, AQ, and STR) fit best, whereas 
the models consisting of old stars (OSP, FR, EXP) fit the poorest. 
Note here that the SSP and AQ models underpredict Mpp/L at 
all cr's. 



As with Fig. [ni the additional tilt of the SSP model 
(i.e. the tilt that is attributable to dark matter and/or 
non-homology, in addition to the stellar populations), is 
small. However, as also noted above, at loga — 1.844, 
the OSP, FR and EXP models overpredict the FP-based 



mass-to-light ratios at the low-cr end by factors 1.45 - 1.7. 
The numerical values are given in Table [S) In contrast, 
the SSP leaves some room (~ 20%) for dark matter. 

6.4. Comparison to dynamical mass-to-light ratios 

A more direct comparison is with dynamical mass-to- 
light ratios in the / and iT-ban d from the Schwarzchil d 
model of the SAURON group (Cappellariitin EOOe'). 
The SAURON results are based on full orbit reconstruc- 
tion using the observed luminosity profile and the ob- 
served line-of-sight velocity distribution. Their method 
thus accounts for spatial and kinematic non-homology 
and so a comparison between a given model and the 
data allow us to determine the fraction of dark matter 
as a function of velocity dispersion. Fig. [13] compares 
the stellar mass-to-l ight ratios of our synth etic galaxies 
with the data from ICappellari et al.l (120061 ) , as a func- 
tion of velocity dispersion^. We find that the slopes of 
the OSP, EXP, and FR models are much flatter than the 
SAURON mass-to-light ratios. The AQ and STR models 
are less flat and the SSP slope is a close match. As noted 
above the masses determined by ICappellari et al.l ()2006l ) 
include both dark and stellar matter, whereas our TMBK 
models predict only the stellar mass-to-light ratio. Note 
tha t aperture effects are sm all: correcting to Re (where 
the ICappellari et al.l ()2006l ) mass-to- light ratios are mea- 
sured) yields changes in log(M,/i) < 0.01. Thus a seri- 
ous problem with the OSP, FR and EXP models is that 
the predict stellar M^/Lj in excess of the total Mdyn/Lj, 
by significant factors (1.93, 1.91 and 1.67, respectively) 
at logcr = 1.844. On this basis, these models are ruled 
out. However, note that even the AQ model has a stellar 
M*/L/ 26% larger than the data. Table |6| outlines the 
offsets between data and model at log(cr) — 1.844. 

Consistent dynamical mass-to-light scaling relations 
are found in other dynamica l stud ies. For example, 
Ivan der Marel fc van DokkumI ( 20071). using a compila - 
tion of data sources including I Cappellari et al.l (|2006D . 
find a steep relation in the B-band, with Mdyn/is oc 
^o.992±o.054 fpj-^jg jg significantly steeper than all models 
except for SSP, for which the slope is 0.773 in the i?-band. 
Their intercept at cr = 100 km/s is log(M/L) = 0.597, 
and is consistent with the SSP model, but, as was the 
case with the /-band, is larger than that of all other 
models at low velocity dispersion. 

From Figure [131 we see that for the SSP case, stel- 
lar populations account for 74% of the total observed 
Mdyn/Li "tilt". Even in the X-band, 63% of the tilt is 
due to stellar population effects. These reduce to 50% 
and 38% for the AQ model. For consistency between our 
models and the SAURON data, we require the baryonic- 
mass to total-mass ratio to scale as /* oc cr""-^^ for SSP, 



and 



for the AQ case. 



The SAURON data are a mixture of cluster and field 
early-types, and it may be argued that it is not well- 
matched to our predictions, which were derived from 
cluster galaxies alone. It is clear from Fig. [131 however, 



^ Note that the ICappellari et aLl II2006I ) data tabulate (Je , the 
velocity dispersion at the effective radius, whereas our velocity 
dispersions are measured at a fixed metric aperture of 0.95 kpc 
radius. For most of the velocity dispersion range. Re is con- 
stant (Fig. [2]l. Since the velocity dispersion aperture correction 
is cr/cTe = (fl^ap/fie)"""®, this yields logcre = logcr — 0.015, which 
is a negligible correction. 
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Fig. 13. — Mass-to-light ratios as a function of velocity dispersion. T he left hand panel show s results in the 7-band; ii'-band results are 
shown on the right. The dynamical mass- to-light ratios from SAURON l lCappellari et an i2006') are shown by the black points (open circles 
are field galaxies, filled circles are cluster galaxies), and the fit to these is indicated by the solid black line. Model fits to the stellar M^/L 
are shown by the colored broken lines (line styles and colors as in Fig. [Sjl. Again, we only fit these lines to log{a) > 1.844, see Section 15.41 
Note that all models except for SSP overpredict the M/L at low a. The models consisting mostly of old stars (OSP, FR, EXP) are the 
poorest fits at low cr. 



TABLE 6 

M/L Offsets at Low Velocity Dispersion 



Model 


Alog(M/L7)=' 


AlogiM/LKT 


AlogiM/Lr)"" 


SSP 


-f0.029 


-0.128 


-0.092 


AQ 


-1-0.099 


-0.066 


-1-0.00 


STR 


-1-0.129 


-0.032 


+0.027 


EXP 


-1-0.222 


-1-0.060 


-f0.157 


FR 


-1-0.281 


-1-0.101 


+0.239 


OSP 


-1-0.286 


-1-0.160 


+0.235 


Offset 


between data and 


model at log((T) 


= 1.844, see Fig. 



I13P Offset between data and model at log{a) = 1.844, see Fig. 

m\ 

that there is no significant difference between cluster and 
field in M^yn/L. At low velocity dispersions (cr < 100 
kms~^) there are 5 galaxies: ICappellari et al.l (|2006f) as- 
sign three of these to clusters (although one of these is 
M32) and the other two to the field. Their H/? values 
range from 2-3 according to their Fig. 16, the upper end 
of which is larger than our median value of 2.05 at the 
same velocity dispersions. Interestingly, one of the high 
H/3, low M/L galaxies is a cluster member. We have 
investigated how the slopes of the SAURON M/L rela- 
tions are affected by dropping between 2 and 4 of the 
SAURON galaxies with the highest H(3 (and hence low- 
est Mdyn/L). We find that this does not significantly 
affect our conclusions: the "mostly-old" class of model, 
when tuned to reproduce the observed spectral indices, 
yields stellar M/L ratios that are too large compared to 
the dynamical values. 

7. EXTENDING THE MODELS 



Our models as described above are necessarily simplis- 
tic. In this section we explore the effects on our results 
of (i) systematics, (ii) altering the IMF, and (iii) adding 
dark matter. 

7.1. Systematic Effects on Mass-to-Light Ratios and 

Colors 

As mentioned in Section r4. 2. 31 the arbitrary choice of 
restricting the star formation of CSFH models to begin 
at a look-back time of istart — 13 Gyr requires investi- 
gation. Again, we choose a look-back time of 10.3 Gyr 
for comparison, which corresponds to z ~ 2. First, we 
note the effect on the colors and mass-to-light ratios of 
the choice of ^start is small. In all models, we find that, 
if the star formation starts at 10.3 Gyr then, at the fixed 
observed values of the line indices, the predicted stellar 
M^,/L values are slightly higher than in the case with 
istart — 13 Gyr, with the effect being largest for the 
bluest bands. This arises because, for example in the 
case of the EXP model, when istart is decreased, then in 
order to fit the observed line indices, t is also reduced 
compared to the istart = 13 Gyr case, with the conse- 
quence that M^/L increases. Other models are affected 
in a similar way. In the r and / bands in the mid-cr range, 
we find a 12% increase in M^/L for the EXP model and 
7% increase for the AQ model. For all models, the slope 
of the Mif/L -a relations are negligibly affected by the 
choice of igtart ■ Colors are also largely unaffected by this 
choice, with the largest offset being for the bluest bands. 
This offset in color towards the blue is < 0.05 in magni- 
tude for all colors given in Tables [5] and [51 for all velocity 
dispersion bins, in all models, compared to the 13 Gyr 
case. This effect can be explained as follows: although 
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^start is now more recent, the constraint on fitting the 
Hne indices forces the fitted 'age' towards older popula- 
tions, as discussed in detail in Section [4.2.3l Because the 
mass-to-light ratios are even higher when star formation 
starts later, this test strengthens the conclusion that the 
'mostly-old' models (EXP, FR, OSP) are a poor fit. 

Next we consider the effect of Balmer line choice, as 
discussed in Section 14.2.21 In that Section, we found 
that, for example, fitting only the H/3 Balmer line had 
the largest effect on the ages^. As with the choice of 
istart, the effect of using only H/3 on the predicted colors 
is negligible. The effect on stellar mass-to-light ratios is 
somewhat greater: for the SSP model, this choice leads 
to a decrease in mass-to-light ratios on the order of 10% 
at log cr ^ 1.85 in the r, / and K bands. For the EXP 
and FR model, however, only a 3% and 1% decreases in 
Mt,/L are found. Thus, this effect is too small to avoid 
our rejection of the FR and EXP models as discussed 
in Section 16. 4( as these models would still overpredict 
the stellar mass-to-light ratios at low velocity dispersion, 
albeit by a slightly lesser amount. 

7.2. Choice of IMF 

In calcula ting stel l ar ma ss-to-light ratios, we have 
adopted the iKroupal (|2001h IMF. Other choices often 
considered in th e literature are the Salpcter (1955) and 
IChabrieJ ()2003f ) IMFs. These models have similar slopes 
at high masses, but differ in the low mass (M < IMq) 
regime: while the Salpeter IMF is a pure power-law, 
the Kroupa and Chabrier IMFs have knees at masses 
of 0.5Mo and IMq respectively. We expect the choice 
of IMF to have little effect on the derived stellar ages 
and metallicities, since most of the light is from main- 
sequence turn-off and red-giant branch stars. For the 
ages considered here (2 Gyr to 13 Gyr), the turn-off 
masses are well above the IMF "knee" . 

The choice of IMF, however, has greater impact on 
the stellar mass-to-light ratios. Using the Salpeter IMF 
would increase the M^,/L ratios by a factor of ^ 1.6 
([Thomas et al.l[2005[ ) and so would increase the conflict 
for all models at low a between the stellar and dynamic 
mass-to-light ratios from the SAURON and FP compar- 
isons. The Chabrier IMF produces marginally lower stel- 
lar mass-to-light ratios compared to a Kroupa IMF, but 
this would be insufficient to reconcile the OSP, EXP and 
FR models with the dynamical mass-to-light ratios. 

More radical IMF rn odifications have been suggested 
by Ivan Dokbm] (|2008D . who proposes an IMF in which 
the critical mass at the knee varies as a function of age, 
with older populations having higher critical masses. It 
is difficult to compare the results of this model in detail, 
as full population synthesi s models have yet to be con- 
structed. Referring to . van DokkumI (|200Sl Fi gure 15, the 
effects on M^/L may be modest, amounting to at most 
a 20% reduction compared to the Chabrier IMF for a 
3 Gyr population. It seems very difficult to reduce the 
high stellar M/L values of the FR and OSP models, as 
most of their population is old, and for older populations 
with a larger critical mass (Mc > O.SA/q), Ad^J L actually 
increases due to remnants. 

^ We note that the other scenarios, such as calibrating the H/3 
hne using all velocity dispersion bins, had a much smaller effect on 
the ages (see Figure|6)l, and hence presumably on the mass-to- light 
ratios, and so we do not discuss those here. 



We conclude that the choice of IMF cannot lower the 
stellar mass-to- light ratios sufficiently to ease the conflict 
with the OSP, FR and EXP models. 

7.3. Effects of Dark Matter and Non-homology 

For simplicity, our models have been assumed to be 
100% steUar, i.e. /» = 1. Various studies have sug- 
gested that, like spiral galaxies, early-type galaxies also 
have non-negligible mass fractions of dark matter (DM) 
in t heir central regions. For example, from lensing stud- 
ies, iBolton et all jlOOl) find /* = 0.62 ± 0.07 within 
the effective radius, for galaxies with 175kms^^ < a < 
400km s~^. It is clear that adding DM will increase the 
total mass-to-light ratios by a factor 1//*. As noted 
above, the OSP, EXP and FR models already over- 
predict the dynamical M/L for low a galaxies, so in- 
cluding DM would only increase the conflict. 

The combined influence of DM and non-homology 
can also be deduced from the tilt of the FP. Follow- 
ing [Boxlan^Kolchi^eOn (|2006D . let us assume a scaling 
Mdyn oc MI^^R^. We define the scaling of the stellar 
mass-to-light as a function of a as {M^,/L) oc cr™, where 
the value of m is given as a function of model and pass- 
band in Tabled The theoretical FP 

Re = {Mdyn/Ly'^a^ihr' (10) 

ZTT 

becomes 

'-(^.)«(^r?^)«<" 

after allowing for stellar populations and dark 
matter/non-homology. If we assume that the dy- 
namical mass-to-light ratio scales with mass M as 
Mdyn/L oc M'^, then the "tilt" T of the FP can be 
expressed as 

T=y^ (12) 

2 + a 

where a is the coefficient of log cr in equation [TlJ We can 
also derive tilts for the cases with only stellar populations 
and no DM (^ = 0, ^ = 0): 



as well as with DM (either ^ or v non-zero). 

Table [7] restates the value of m for each model from 
Table [5] in the r-band. We find that all of the models 
with the exception of the SSP are too steep (too high 
a, close to the virial relation) in comparison to the ob- 
served FP. We can then adjust the DM/non- homology 
scaling (the value of fi or v) so as to match the observed 
FP slope in Figure [TT] This yields the values in the sec- 
ond and third columns of Table [T] For the SSP model, 
the DM/non-homology scaling is essentially zero. For 
the other models, stronger scaling is required: for AQ 

= 0.06, or V — 0.150. As a result, we find that in the 
SSP case, essentially all of the r-band FP tilt is due to 
stellar populations. For the AQ model, 2/3 of the tilt is 
due to stellar populations and 1/3 is due to DM and/or 
non-homology. Note that because of the age-dependence 
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TABLE 7 

Effect of Dark Matter on t-band FP Tilt 



Model 














SSP 


0.662 


0.001 


0.003 


0.198 


0.200 


0.993 


AQ 


0.455 


0.065 


0.150 


0.128 


0.197 


0.652 


STR 


0.391 


0.085 


0.190 


0.108 


0.194 


0.558 


EXP 


0.235 


0.140 


0.310 


0.062 


0.194 


0.322 


FR 


0.174 


0.165 


0.355 


0.045 


0.194 


0.234 


OSP 


0.196 


0.155 


0.340 


0.052 


0.194 


0.270 



" (M*/L)r oc cr™. From table |9f required in eq. [TTI fwith 
V = 0) to reproduce the log(o') coefficient in the observed Coma 
FP.'^ same as (b), but with fi = 0'^ M^/L oc W^' : tilt due to stellar 
populations'^ M/L oc M^: tilt of the FP* Tt/T is the fraction of 
the tilt due to stellar populations 



of the K-haiid stellar M/L, for the SSP case we find that 
most of the predicted /C-band tilt is due to stellar pop- 
ulation effects. In the case of the AQ model, we expect 
that half of the iiT-band tilt is due to stellar population 
effects. 

If we repeat this exercise using the dynamical Mp p / L 
values calculated from the FP and shown in Fig. [T^l we 
can determine the tilt assuming that dark matter and 
non-homology depend directly on a. For the SSP case, 
we find that /* oc ct""'^'^ and for the AQ case the expo- 
nent is —0.34. To extend these scaling relations to the 
case of non-homology, replace with the combination 
/*(c2/5). These scalings are slightly flatter than the cor- 
responding values derived from the SAURON compari- 
son: oc (T~°'^^ and oc (T~°-^^, respectively. Since the 
SAURON results are free of the effects of non-homology, 
the comparison between the FP and SAURON scalings 
then yields C2 oc cr^"^ for both SSP and AQ cases. This 
non-homology is weak, and is not a major driver of the 
tilt of the FP. 

8. DISCUSSION 

8.1. Summary of Model Fits 

We consider each model separately, and summarize the 
comparison with the observables. The OSP case is al- 
ready ruled out by the Balmer linestrengths, but we note 
that it is also a poor fit to the colors, to the FP and dy- 
namical M/L. 

FR: The 2% frosting model consists mostly of old stars. 
While this model is a good fit to the slope of the 
u — g color-cr relation, it is a poor fit to the slopes of 
the g — i and particularly the i — K color-cr relation. 
It overpredicts the FP-based r-band mass-to-light 
ratios, and those from SAURON by significant fac- 
tors (1.7 - 1.9). We reject this model. 

EXP: This model also consists mostly of old stars, with 
the young ages arising from the exponential tail of 
young stars. This model, like FR, is a good fit to 
the slope of the u — g diagram, it is too flat flt 
in g — i. Like the FR model, the stellar M/L also 
overpredicts the FP-based dynamical r-band mass- 
to-light ratios, and those from SAURON by factors 
1.45 & 1.65, respectively. Given that we expect 
there to be some dark matter present in RSGs, we 
rank this model as "disfavored" . 



AQ: The AQ model consists of old through intermedi- 
ate age stars. In the color diagrams, its fit is better 
than EXP but poorer than SSP. In terms of mass- 
to-light ratio, it marginally overpredicts the FP and 
dynamical mass-to-light ratios, by factors 1.09 and 
1.25. Given the uncertainties, we rank this model 
as "acceptable", although we note that the impli- 
cation of this model is that there is little DM at 
low velocity dispersions. 

STR The behavior of this model is similar to the AQ 
model in most respects. 

SSP: This model is the best fit to the linestrengths. It 
also fits the slopes of the color diagrams the best 
overall: it is somewhat too steep in u-g, but the 
expectation that there may be small amounts of 
dust affecting the observed u-g colors at low a goes 
in the right sense to account for this difference. It 
is the best flt in 5 — i, and aside from the rejected 
OSP model has the closest fit irii — K. Its predicted 
stellar mass-to-light ratios are always comfortably 
below the FP-based by factors ^ 0.8 and appear 
to be similar to the dynamical mass-to-light ratios 
measured by the SAURON team. We rank this 
model as "good" . 

In summary, the models that fit best to the faint end 
of the FP, and to the SAURON mass-to-light ratios, are 
the SSP and AQ models. The models which are "mostly- 
old" (EXP, FR, OSP) faff to fit the faint end of the FP 
or the low-(T end of the SAURON M/L ratios by factors 
of - 1.45-1.7. 

8.2. Comparison of the predicted and observed scatter 

For all of the models, with the exception of SSP, 
the predicted scatter from stellar population effects is 
well within the scatters observed in the color-a, M/L-cr, 
and FP relations. The SSP model, however, overpre- 
dicts the scatter in the M/L and, to a lesser extent, in 
the FP. This is a direct result of assuming that age is 
strongly anti-correlated with metallicity at a given ve- 
locity dispersion (Section 14. 1.2p . Recall that the scatter 
in stellar population parameters is constrained by the 
scatter in linestrengths. When age and metallicity are 
anti-correlated as we have assumed, their effects on the 
predicted linestrength tend to partially cancel. Conse- 
quently, the observed scatter in a given linestrength in- 
dex can be fit with a larger scatter in age and a larger 
anti-correlated scatter in metallicity than would be the 
case if these two parameters were assumed uncorrelated. 
The same "conspiracy" keeps the scatter in colors low. 
Mass-to-light ratios are more sensitive to age, however, 
and so as a result, the scatter in M^/L is larger than it 
would be if we had assumed that age and metallicity were 
uncorrelated when fitting the linestrength scatters. As a 
result, this does not invalidate the SSP model itself, since 
it may be possible to fit the line strengths adequately 
with a weaker age-metallicity anti-correlation^. For the 
other CSFH models, the stellar M/L ratios are largely 
dominated by older stars and so have considerably less 
scatter. 

^ See ISmith et al.l l)2009al ) for an alternative approach to con- 
straining the anti-correlation between age and metallicity. 
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8.3. Timescales of Star Formation Histories of 
Red-Sequence Galaxies 

The results here argue that star formation histories of 
RSGs were either a short SSP-like burst or of somewhat 
more extended duration as in the AQ model. It is difficult 
to distinguish between these two cases based on this low- 
redshift data alone. 

One method, although it is somewhat uncertain, is via 
a-enhancement ratios. iThomas et al.l (|2005ll have ar- 
gued that, in order to obtain high a-enhancement ra- 
tios, the duration of star formation must be less than 
^ 1 Gyr. Specifically, they propose a relation [a/Fe] = 
1/5 - l/61og(At), where At is the fuU-width-at-half- 
maximum of a Gaussian star formation history. For 
RSGs of all masses, [a/Fe] > 0.2 (Fig. Hand El) so this 
implies At < 1 Gyr. Thus based on a-enhancement ra- 
tios and the lThomas et all (|2005D models, the SSP would 
be preferred over the AQ model, for which the duration 
of star formation can extend over a range of up to ~ 10 
Gyr at cr ~ 70 kms~^. However, it must be emphasized 
that the translation between [a/Fe] and star- formation 
timescale depends on many assumptions about the chem- 
ical enrichment histories of early- type galaxies, that re- 
main uncertain at present. 

Of course, more complex star formation histories than 
the ones described here are possible. For example, a plau- 
sible model is one in which the bulge component formed 
in a short SSP-like burst, whereas the disk component 
was quenched. Clearly the results from this model would 
lie somewhere between the SSP and AQ cases. 

8.4. Implications for the tilt of the FP 

We have shown that the best-fitting model is the SSP, 
and that, as a result of the scaling of SSP age with ve- 
locity dispersion, there is a strong trend of stellar mass- 
to-light ratio with velocity dispersion. Consequently, we 
find that 75%-100% of the r//-band FP tilt is due to 
stellar population variations as a function of a. 

A strong SSP-age-driven scaling has generally not been 
assumed in previous analyses of the FP, leading to the 
conclusion that the "extra" tilt, e.g. due to DM varia- 
tions, is much stronger than we find here. For example, 
iTrujillo et al.l (|2004l ). assumed an exponential star for- 
mation history model, and concluded that only 1/4 the 
tilt was due to stellar populations. In fact, we would 
have come to a similar value (30%), had we adopted the 
EXP model (see Table [71). The EXP model is, however, 
inconsistent with the observed dynamical mass-to-light 
r atios, as we have shown (Sec tion 16. 4p . 

iPadmanabhan et al.l ()2004f ) based t heir scaling anal- 
yses o n the star formation models of iKauffman n et al] 
(l200l . which are combinations of the EXP and FR 
models, and found a strong scaling of dynamical mass 
with stellar rn a ss. A similar conclusion was reached by 
iGallazzi et all (l200l who fo und Mdy„/M* cx M^-^s, 
and iHvde fc Bernardil (|2008H) . who found Mdy„/M* cx 
-^dyr/i both ba sed on models similar to those of 
IKauffmann et al] (|2003l ). As a result the assumed star 
formation history, these studies will overestimate the 
stellar mass at low masses and hence will overestimate 
the strength of the DM scaling. For the SSP model, there 
is essentially no FP tilt due to DM. For the AQ model, 
the scaling is weak: M^y^/M^ cx Af^ °^. 



These results have important consequences for esti- 
mating the contributions of mergers and dissipation in 
forming the FP. For exampl e, dry mergers predict a scal- 
ingMHyn/M, cx M°-i^- ° '^5 (iBovlan-Kolchin et al.l2 0061. 
iRobertson et al.l ()2006l ) studied dissipational mergers ar- 
gued that the tilt of the FP is due to a varying amount of 
dissipation in mergers along the sequence as a function 
of mass. If this interpretation is correct, a consequence 
of our study is that there is little, if any, variation in 
dissipation as a function of mass. 

One might assume that the tilt of the FP in the near- 
infrared would be less sensitive to the effects of stellar 
populations, and hence that the ii'-band FP tilt should 
reflect primarily the effects of dark matter and non- 
homology. However, this is approximately valid only 
in mostly-old SFH models, which cannot simultaneously 
match line indices and the absolute mass-to-light ratio. 
For SSP or AQ models, the K-band tilt can be strongly 
affected by stellar population variations. 

8.5. Implications for stellar masses and the stellar mass 
density of the Universe in red galaxies 
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Fig. 14. — M/Lk vs Mk- Best fit lines are calculated for 
those galaxies below Mk < —21, where the data are complete is 
log(cr). We have overlaid the contours for the SSP and EXP case, 
to illustrate that the scatter in \og{cr) is diagonal, and also to show 
the very different scatters in log(Af/Lif) between models. The 
quoted zero-points (z.p.) are measured at Mk = —24 (L*). Line 
styles and colors as in Fig. |3] Green and blue contours are for the 
SSP and EXP populations, and levels are as in Fig. |4] Notice the 
larger scatter in the SSP case, which decreases the mean mass-to- 
light ratio. 

We have found that models that are "mostly-old" are 
ruled out: for small mass systems, their stellar M/L ra- 
tios exceed the total dynamical M/L ratios by factors 
of ^ 1.45 — 1.7. This class of models, which includes 
exponential star formation histories and "frosting" sce- 
narios, have often be en used to model t he star formation 
histories of galaxies (|Cole et al.l 120011 : IKauffmann et all 
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[200I iBell et alJ[200l iGallazzi et al\\20M) . While this 
may be a reasonable model for late- type or blue galaxies, 
we have found that it is not a good fit for RSGs. Con- 
sequently, the stellar masses of red galaxies have been 
overestimated by this assumption. 

If we assume a different star formation history model, 
we will arrive at different conclusions, even at bright 
magnitudes. This is because the dependence of M^/Lk 
on age is quite strong for e.g. the SSP model, and hence 
the scatter in M^:/Lk is also quite large. Thus many 
galaxies which are bright in K are actually galaxies with 
relatively small stellar mass that are quite young, boost- 
ing their if -band luminosity and lowering their mass-to- 
light ratio. Figure [T3] shows the if -band mass-to-light 
ratio for the different star formation models. The ratios 
of the mean if -band mass-to- light ratios at Mj^ = — 24 
(the characteristic A'-band luminosity) are 0.73:0.82:1 for 
SSP:AQ:EXP. Thus for the SSP case, the luminosity den- 
sity in the early-type galaxies or RSGs as calculated from 
the if-band should be reduced by a factor of ^ 0.73. 
This woul d redu ce the stellar mass density calculated by 
IBell et al.l (j2003D from 3.2 in units of IQ-^/iMq Mpc'^ 
to 2.3 in the same units, and would reduce the fraction of 
stellar mass in early types (assuming that late-types are 
unaffected) from 60% to 50%. Repeating this calculation 
in the g-band, for example, the reduction in stellar mass 
in RSGs is similar to that found for morphologically- 
selected early-types in the if -band: a factor ~ 0.7. Thus 
the total fraction of all stell ar mass wh i ch res ides on the 
red sequence, as derived bv IBell et al.l (|2003f ). would be 
reduced from 70% to 60%. 

9. CONCLUSIONS 

We have constructed 6 parametric star formation mod- 
els, each representing a distinct star formation his- 
tory. The modeled Lick indices were fit to observed 
linestrength data in the Shapley supercluster, and sets 
of best-fit model parameters were constructed for each 
velocity dispersion bin for each model. Synthetic clus- 
ter populations were then constructed based on the sta- 
tistical distributions of each stellar population parame- 
ter, and mass-to-light ratios were derived from the stellar 
populations. 

Our main results are as follows: 

1. From the fits to linestrength indices, we find that 
in all models there exists a "downsizing" trend in 
luminosity weighted-age age for galaxies with a > 
70km s^^. Lower-tT galaxies also tend to be less 
metal rich and a-enhanced. We find that, based on 
the fits to the Lick indices alone, the OSP model is 
strongly rejected. This agrees with previous studies 
(Nclan'ct al. 2005). The EXP and FR models are 
poorer fits to the linestrength data than the SSP 
model, but cannot rejected from this test alone. 

2. The linestrength data suggest that the "downsiz- 
ing" effect stops at tr ~ 70 kms~^ below which 



luminosity-weighted stellar ages either stop de- 
creasing or begin to increase. This is not a selection 
effect, which would tend to bias the data in the op- 
posite sense. 

3. The derived stellar mass-to-light ratios from the 
models were compared to dynami cal mass-to-light 
ratios from the FP and from ICappellari et al.l 
(|2006f ). The "mostly-old" class of star formation 
models (OSP, EXP and FR) have Kroupa-IMF stel- 
lar mass-to-light ratios that are larger than the 
observed dynamical mass-to-light ratios for low-cr 
galaxies, by factors 1.45-1.7. For these models, the 
addition of dark matter would make the agreement 
worse, thus these "mostly-old" models are ruled 
out. We conclude that the ^ 6 Gyr luminosity- 
weighted ages of low-mass RSGs do indeed reflect 
the "intermediate" ages of the bulk of the stellar 
population, as parametrized by the SSP or AQ 
mdoels. A hybrid of the SSP and AQ, in which, 
for example, the bulge of the system formed in a 
short intense burst, while the disk was quenched at 
a later time, would likely fit all of the data well. 

4. We have shown that, for the "intermediate-age" 
scenarios, the scaling of stellar mass-to-light ratio 
with mass is strong, both in the optical bands and 
in the if -band. Consequently, stellar populations 
explain most of the FP tilt in the optical bands, 
and, for the SSP case, at least 50% of the tilt in 
the if -band. This leaves much less room for tilt 
due to variations in the dark-to-stellar mass along 
the red sequence. 

5. For the SSP model, the stellar masses are consid- 
erably reduced, particularly at low velocity disper- 
sions, compared to the star formation histories as- 
sumed in previous studies. As a consequence, the 
stellar mass density in red galaxies may be as much 
as 30% lower than previously assumed. 

A further test of the star formation histories of RSGs 
is via the redshift evolution of their red and blue frac- 
tions, of their dwarf-to-giant ratios and of their color- 
magnitude relation slopes and zero-points. In particular, 
the AQ and SSP models will show strong evolution from 
the blue cloud to the red-sequence, whereas the "mostly- 
old" models will have a well-populated red-sequence at 
intermediate lookback times (z ~ 0.5). We will explore 
these predictions in future work. 

We thank Claudia Maraston for providing tables of 
mass-to-light ratios based on empirical stellar spectra 
in advance of publication. MJH acknowledges support 
from NSERC. RJS is supported by STFC rolling grant 
PP/C501568/1 'Extragalactic Astronomy and Cosmol- 
ogy at Durham 2005-2010'. 



APPENDIX 

TABLES OF SCALING RELATIONS 

Tables ISI and |5| summarize the slopes and intercepts as a function of velocity dispersion for Johnson-Cousins, infrared 
and SDSS passbands for each star formation history model. Zero points are measured at log(cr) — 2.00. 
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Mass-To-Light Ratio, Magnitude, color 



TABLE 8 

Velocity Dispersion Relations By Model (Johnson-Cousins 



Band (X) 


^ _ dlog(M/Lx)a. 


z.p.t* 


r.m.s. 


dMx a 


z.p.t* 


r.m.s. 


color 


d(X-X')a 




r.m.s. 


d log(o-) 


dlog(CT) 


dlog(CT) 






Single Stcll 


ar Populat 


ion Model (SSP) 










U 


0.956 


0.757 


0.135 


-3.666 


-18.283 


0.661 


U-B 


0.457 


0.571 


0.041 


B 


0.773 


0.581 


0.125 


-4.123 


-18.855 


0.649 


B-V 


0.199 


0.980 


0.019 


V 


0.693 


0.449 


0.120 


-4.322 


-19.835 


0.643 


V-R 


0.099 


0.475 


0.009 


R 


0.654 


0.423 


0.118 


-4.421 


-20.310 


0.640 


R-I 


0.093 


0.613 


0.010 


I 


0.617 


0.313 


0.116 


-4.513 


-20.923 


0.638 


I-J 


0.231 


1.107 


0.046 


J 


0.524 


0.047 


0.128 


-4.744 


-22.030 


0.653 


J-H 


0.054 


0.674 


0.008 


H 


0.503 


-0.095 


0.128 


-4.798 


-22.704 


0.653 


H-K 


0.052 


0.224 


0.017 


K 


0.482 


-0.168 


0.133 


-4.851 


-22.928 


0.660 










Abruptly Quenched Model (AQ) 


U 


0.726 


0.799 


0.082 


-4.231 


-18.191 


0.606 


U-B 


0.422 


0.574 


0.052 


B 


0.557 


0.621 


0.069 


-4.653 


-18.764 


0.597 


B-V 


0.186 


0.981 


0.022 


V 


0.483 


0.489 


0.064 


-4.839 


-19.745 


0.594 


V-R 


0.087 


0.477 


0.012 


R 


0.448 


0.462 


0.062 


-4.926 


-20.222 


0.593 


R-I 


0.076 


0.615 


0.012 


I 


0.418 


0.352 


0.061 


-5.002 


-20.837 


0.592 


I-J 


0.227 


1.114 


0.060 


J 


0.327 


0.082 


0.072 


-5.230 


-21.952 


0.601 


J-H 


0.041 


0.677 


0.011 


H 


0.311 


-0.061 


0.074 


-5.270 


-22.628 


0.602 


H-K 


0.048 


0.227 


0.021 


K 


0.292 


-0.135 


0.080 


-5.318 


-22.856 


0.607 
















Strangulation M 


odel (STR) 












U 


0.664 


0.819 


0.062 


-4.353 


-18.132 


0.596 


U-B 


0.421 


0.573 


0.048 


B 


0.496 


0.642 


0.049 


-4.775 


-18.705 


0.590 


B-V 


0.191 


0.979 


0.020 


V 


0.419 


0.510 


0.044 


-4.966 


-19.684 


0.587 


V-R 


0.086 


0.477 


0.010 


R 


0.385 


0.483 


0.042 


-5.052 


-20.161 


0.587 


R-I 


0.073 


0.616 


0.011 


I 


0.356 


0.373 


0.042 


-5.125 


-20.777 


0.587 


1-.) 


0.235 


1.111 


0.057 


J 


0.262 


0.105 


0.052 


-5.360 


-21.888 


0.594 


J-H 


0.039 


0.677 


0.009 


H 


0.246 


-0.038 


0.054 


-5.399 


-22.564 


0.595 


H-K 


0.051 


0.226 


0.018 


K 


0.226 


-0.112 


0.060 


-5.450 


-22.790 


0.599 











Exponential Stcir Formation Hate Model (EXP) 



u 


0.478 


0.896 


0.063 


-4.694 


-17.931 


0.597 


U-B 


0.380 


0.583 


0.063 


B 


0.326 


0.715 


0.040 


-5.073 


-18.513 


0.582 


B-V 


0.168 


0.982 


0.024 


V 


0.259 


0.582 


0.031 


-5.241 


-19.496 


0.578 


V-R 


0.072 


0.481 


0.013 


R 


0.230 


0.554 


0.027 


-5.313 


-19.977 


0.577 


R-I 


0.063 


0.621 


0.013 


I 


0.205 


0.442 


0.025 


-5.376 


-20.597 


0.575 


I-J 


0.230 


1.107 


0.071 


J 


0.113 


0.175 


0.032 


-5.607 


-21.705 


0.577 


J-H 


0.040 


0.679 


0.013 


H 


0.096 


0.031 


0.036 


-5.647 


-22.384 


0.578 


H-K 


0.048 


0.225 


0.020 


K 


0.077 


-0.043 


0.042 


-5.694 


-22.609 


0.581 











Frosting Model (FR)~ 



u 


0.395 


0.972 


0.069 


-5.123 


-17.743 


0.596 


U-B 


0.341 


0.613 


0.062 


B 


0.259 


0.779 


0.044 


-5.464 


-18.357 


0.581 


B-V 


0.159 


0.996 


0.025 


V 


0.195 


0.640 


0.034 


-5.623 


-19.353 


0.577 


V-R 


0.062 


0.489 


0.012 


R 


0.171 


0.609 


0.030 


-5.684 


-19.842 


0.576 


R-I 


0.048 


0.630 


0.014 


I 


0.151 


0.493 


0.025 


-5.733 


-20.472 


0.574 


I-J 


0.201 


1.123 


0.063 


J 


0.071 


0.220 


0.011 


-5.934 


-21.594 


0.572 


J-H 


0.014 


0.689 


0.010 


H 


0.066 


0.072 


0.013 


-5.947 


-22.283 


0.573 


H-K 


0.026 


0.233 


0.017 


K 


0.055 


-0.005 


0.019 


-5.973 


-22.517 


0.574 











Old Single Stellar Population Model (OSP) 



u 


0.451 


0.960 


0.073 


-4.903 


-17.774 


0.603 


U-B 


0.423 


0.584 


0.069 


B 


0.282 


0.779 


0.046 


-5.326 


-18.358 


0.586 


B-V 


0.159 


0.993 


0.026 


V 


0.218 


0.641 


0.035 


-5.486 


-19.352 


0.582 


V-R 


0.068 


0.484 


0.011 


R 


0.191 


0.612 


0.031 


-5.553 


-19.835 


0.580 


R-I 


0.069 


0.620 


0.011 


I 


0.164 


0.500 


0.027 


-5.622 


-20.456 


0.579 


I-J 


0.303 


1.057 


0.048 


J 


0.043 


0.253 


0.008 


-5.926 


-21.512 


0.577 


J-H 


0.059 


0.672 


0.011 


H 


0.019 


0.112 


0.004 


-5.985 


-22.185 


0.576 


H-K 


0.062 


0.210 


0.010 


K 


-0.006 


0.044 


0.004 


-6.047 


-22.394 


0.577 











Measured with galaxies with log(iT) > 1.844*^ Measured at log{a) = 2.00 
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TABLE 9 

Mass-To-Light Ratio, Magnitude, color - Velocity Dispersion Relations By Model (SDSS [AB]) 



R-inrl ^\^ 1 m — l°s{*f/ix ) a ™ „ b r tti - dMx a r, „ h r m - r-r,1nr- <i(X-X')i^ „ b r m - 

tsana m dlog(CT) ^-P- dlog(<T) ^-P- ^^^^^^ dlog(CT) ^-P- r.m.s. 




Single Stellar Population Model (SSP) 


u 

g 

r 
i 

z 


0.958 0.735 0.136 
0.740 0.528 0.123 
0.662 0.373 0.119 
0.628 0.289 0.117 
0.585 0.191 0.118 


-3.529 -17.553 0.667 
-4.075 -19.360 0.651 
-4.270 -20.198 0.646 
-4.354 -20.518 0.643 
-4.463 -20.814 0.645 


u-g 
g_r 
r-i 
i-z 
z-J 


0.545 1.807 0.050 
0.196 0.8.38 0.018 
0.084 0.320 0.009 
0.109 0.296 0.015 
0.151 1.220 0.034 


Abruptl 


y Quenched Model (AQ) 


u 

g 
r 
i 
z 


0.727 0.777 0.082 
0.527 0.568 0.067 
0.455 0.412 0.062 
0.428 0.328 0.061 
0.388 0.228 0.063 


-4.151 -17.456 0.608 
-4.653 -19.266 0.597 
-4.833 -20.106 0.594 
-4.900 -20.428 0.594 
-5.000 -20.726 0.595 


u-g 
g-r 
r-i 
i-z 
z-J 


0.502 1.811 0.061 
0.179 0.840 0.022 
0.068 0.322 0.010 
0.100 0.298 0.021 
0.153 1.226 0.044 


Strangulation Model (STR) 


u 

g 
r 
i 

z 


0.664 0.797 0.061 
0.464 0.589 0.047 
0.391 0.434 0.043 
0.366 0.349 0.042 
0.325 0.250 0.043 


-4.253 -17.403 0.586 
-4.753 -19.212 0.579 

-4.934 -20.051 0.577 
-4.999 -20.374 0.577 
-5.101 -20.671 0.578 


u-g 

g-r 

r-i 
i-z 
z-J 


0.500 1.809 0.057 
0.181 0.839 0.020 

0.065 0.323 0.009 
0.102 0.298 0.020 
0.158 1.223 0.041 


Exponential Star Formation Rate Model (EXP) 


u 

g 
r 
i 

z 


0.477 0.876 0.064 
0.297 0.663 0.036 
0.235 0.505 0.028 
0.214 0.418 0.025 
0.174 0.319 0.023 


-4.830 -17.201 0.589 
-5.280 -19.023 0.573 
-5.435 -19.868 0.570 
-5.488 -20.196 0.570 
-5.588 -20.494 0.569 


u-g 

g-r 
r-i 
i-z 
z-J 


0.450 1.822 0.075 
0.155 0.844 0.024 
0.053 0.328 0.010 
0.100 0.298 0.026 
0.154 1.220 0.050 


Frosting Model (FR) 


u 

g 
r 

i 
z 


0.393 0.952 0.069 
0.231 0.724 0.040 
0.174 0.561 0.031 
0.159 0.470 0.027 
0.124 0.367 0.018 


-5.090 -17.003 0.593 
-5.494 -18.863 0.578 
-5.637 -19.721 0.574 
-5.675 -20.058 0.573 
-5.761 -20.365 0.572 


u-g 
g-r 
r-i 
i-z 
z-J 


0.404 1.860 0.073 
0.143 0.859 0.024 
0.0.38 0.337 0.011 
0.086 0.307 0.025 
0.1.34 1.229 0.044 


Old Single Stellar Population Model (OSP) 


u 

g 
r 

i 

z 


0.454 0.939 0.074 
0.255 0.724 0.041 
0.196 0.563 0.032 

0.176 0.175 0.029 
0.12 i 0.382 0.020 


-4.873 -17.033 0.599 
-5.370 -18.859 0.581 
-5.519 -19.710 0.578 

-5.569 -20.011 0.577 
-5.699 -20.321 0.57(5 


u-g 
g-r 

r-i 

i-z 
z-,1 


0.498 1.826 0.082 
0.149 0.851 0.024 

0.050 0.330 0.008 

0.129 0.284 0.021 
0.201 1.182 0.032 


^ Measured with galaxies with log(iT) > 1.844'' Measured at log(iT) = 2.00 
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